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ABSTRACT 

• We present a new evolutionary stellar population synthesis model, which pre- 

dicts SEDs for single-age single-metallicity stellar populations, SSPs, at resolution 
1.5A (FWHM) in the spectral region of the near-IR Call triplet feature. The main 
ingredient of the model is a new extensive empirical stellar spectral library that has 
been recently presented in Cenarro et al. (2001a, b; 2002), which is composed of more 
than 600 stars with an unprecedented coverage of the stellar atmospheric parameters. 

Two main products of interest for stellar population analysis are presented. The 
first is a spectral library for SSPs with metallicities — 1.7<[Fe/H]<-|-0.2, a large range 
of ages (0.1-18 Gyr) and IMF types. They are well suited to model galaxy data, since 
the SSP spectra, with flux-calibrated response curves, can be smoothed to the resolu- 
, tion of the observational data, taking into account the internal velocity dispersion of 

the galaxy, allowing the user to analyze the observed spectrum in its own system. We 
also produce integrated absorption line indices (namely CaT*, CaT and PaT) for the 
same SSPs in the form of equivalent widths. 

We find the following behaviour for the Call triplet feature in old-aged SSPs: i) 
|— ) \ the strength of the CaT* index does not change much with time for all metallicities 

for ages larger than ^3 Gyr, ii) this index shows a strong dependence with metallicity 
for values below [M/H] ~ —0.5 and in) for larger metallicities this feature does not 
show a significant dependence either on age or on the metallicity, being more sensitive 
to changes in the slope of power-like IMF shapes. 
' The SSP spectra have been calibrated with measurements for globular clusters of 

Armandroff & Zinn (1988), which are well reproduced, probing the validity of using 
the integrated Can triplet feature for determining the metallicities of these systems. 
Fitting the models to two early-type galaxies of different luminosities (NGC 4478 and 
NGC 4365), we find that the Can triplet measurements cannot be fitted unless a very 
dwarf-dominated IMF is imposed, or if the Ca abundance is even lower than the Fe 
abundance. More details can be found in Cenarro et al. (2003). 

Key words: stars: fundamental parameters — globular clusters: general — galaxies: 
abundances — galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: 
stellar content 
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1 INTRODUCTION 

This is the fourth and last paper in a series devoted to un- 
derstand the behavior of the near-infrared Can triplet fea- 
ture in stars and in stellar populations. The ultimate aim of 
this work is to use the strength of the Call triplet lines in 

* E-mail: vazdekis@ll.iac.es 



this spectral range to investigate the stellar content of early- 
type galaxies, but the results are sufficiently general to be 
used in other areas (starburst and active galaxies, globular 
clusters, or stellar astrophysics) as well. In Cenarro et al. 
(2001a) (hereafter Paper I) we presented a new empirical 
stellar spectral library, mostly observed at the 1 m Jacobus 
Kapteyn Telescope (JKT) at the Observatorio del Roque de 
Los Muchachos, La Palma. The library is composed of 706 
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stars covering the Can triplet feature in the near-infrared. 
In that paper we provided a new index definition for this fea- 
ture (CaT*) that, among other advantages, minimizes the 
effects of the Paschen series. For the stars in this library 
we derived a set of homogeneous atmospheric parameters in 
Cenarro et al. (2001b) (hereafter Paper II). In Cenarro et 
al. (2002) (hereafter Paper III) we described the behaviour 
of the Can triplet through empirical fitting functions, which 
relate the strength of this feature to the stellar atmospheric 
parameters. In the current paper, we make use of these func- 
tions and the spectra of the stellar library to predict both 
the strength of the Call triplet feature and spectral energy 
distributions, SEDs, in the range AA8348.85-8950.65A for 
single-age, single-metallicity stellar populations (SSPs) by 
means of evolutionary stellar population synthesis model- 
ing. 

Traditionally, elliptical galaxies have been thought to 
be a uniform class of objects, with global properties chang- 
ing smoothly with mass and hosting old and coeval stellar 
population. However, over the last decade, a growing body 
of evidence is indicating that the formation processes and 
star formation histories of, at least, an important fraction 
of early- type galaxies are more complex and heterogeneous. 
The apparent age spread among elliptical galaxies (Gonzalez 
1993; Faber et al. 1995; j0rgensen 1999), the distribution of 
element abundances (Worthey 1998; Peletier 1999; Trager et 
al. 2000a) and the interpretation of the scaling relations (like 
the colour-magnitude or Mg2-(J relations (Bower, Lucey & 
Ellis 1992; Bender, Burstein & Faber 1993; Colless et al. 
1999; Terlevich et al. 1999; Kuntschner 2000, hereafter K00; 
Trager et al. 2000b, hereafter TOO; Proctor & Sansom 2002, 
hereafter PS02), are some of the main issues in the present 
debate about the evolutionary status of early-type galaxies. 
These studies have been possible thanks to the comparison 
of observed data to the predictions of the, so called, stel- 
lar population synthesis models. These models make use of 
a theoretical isochrone, or H-R diagram, convert isochrone 
parameters to observed ones, assuming empirical or theoret- 
ical prescriptions, and finally integrate along the isochrone 
assuming an initial mass function, IMF (e.g. Tinsley 1980; 
Bruzual 1983; Arimoto & Yoshii 1986). 

A large part of the discussion presented above is based 
on observations in the Lick/IDS system (Worthey et al. 
1994, hereafter W94, and references therein), a set of 21 
absorption line indices from 4100 to 6300 A. Based on high 
quality data in that wavelength region, and stable stellar 
population models we are left with the puzzling result that 
there is large scatter in the luminosity-weighted mean age 
of the elliptical galaxies (e.g. TOO). On the other hand the 
colour-magnitude relation at z = 0.8 continues to exist (El- 
lis et al. 1997; Stanford et al. 1997), indicating that at least 
a large fraction of elliptical galaxies in the local universe 
are coeval. However an unambiguous analysis of the stellar 
populations has been hampered by the fundamental age- 
metallicity degeneracy, i.e. the two effects cannot be fully 
separated in the integrated spectrum of a composite stel- 
lar population. Some newer, blue indices (Jones & Worthey 
1995; Worthey & Ottaviani 1997, hereafter W097; Vazdekis 
& Arimoto 1999, hereafter VA99; Gorgas et al. 1999) could 
contribute to alleviate this situation. 

In this set of papers we have increased the wavelength 
region that can be studied, adding the tools to analyze a 



few important absorption lines in the near-infrared. In fact, 
by extending the spectral coverage, we increase our stellar 
population analysis constraining power, since the contribu- 
tion of different types of stars to the total luminosity varies 
as a function of the spectral range. Furthermore, contrary 
to the ultraviolet region, where a large fraction of the light 
can originate from just a few stars (e.g. Ponder et al. 1998), 
most of the stars contribute to the near-infrared Ca II triplet 
lines. Thanks to the large stellar library and state-of-the-art 
evolutionary models provided in this series of papers, ob- 
servers will now be able to analyze more accurately their 
galaxy Ca II triplet measurements, and compare them with 
the models that fit in the blue. 

The Call triplet is one of the most prominent features 
in the near-IR spectrum of cool stars and its potential to 
study the properties of stellar populations has been exten- 
sively acknowledged in the literature. For instance, the Call 
triplet strength has been found to correlate with globular 
cluster metallicity and, therefore, it has been proposed as a 
metallicity indicator for old and coeval stellar populations 
for the metallicity regime typical of galactic globular clus- 
ters (Armandroff & Zinn 1988, hereafter AZ88). This rela- 
tion has been recalibrated by Rutledge, Hesser & Stetson 
(1997), which review the different methods to measure clus- 
ter metallicities using this feature. Terlevich, Diaz & Ter- 
levich (1990a) found that active galaxies exhibit Call triplet 
strengths equal or larger than those found in normal el- 
lipticals, which they interpreted as due to the presence of 
red supergiant stars in the central regions of these galax- 
ies. The same approach has been followed by a number of 
authors (Forbes, Boisson & Ward 1992; Garcia- Vargas et 
al. 1993; Gonzalez Delgado & Perez 1996ab; Hcckman et al. 
1997; Perez et al. 2000). The suggested gravity sensitivity 
of the Can triplet has also been proposed to constrain the 
dwarf/giant ratio in early-type galaxies (e.g. Cohen 1978; 
Faber & French 1980; Carter, Visvanathan & Pickles 1986; 
Alloin & Bica 1989). This feature has been measured in ex- 
tensive galaxy samples (e.g. Cohen 1979; Bica & Alloin 1987; 
Terlevich et al. 1990b; Houdashelt 1995). Among the most 
interesting results derived from these studies is the fact that 
the Ca II triplet strength does not seem to vary much among 
early- type galaxies of different types, colours and luminosi- 
ties. This result is at odds with the strong metallicity corre- 
lation found for the globular clusters. We refer the reader to 
§ 2 of Paper I for a review of previous works on the subject. 

A reliable analysis of the Can triplet measurements in 
integrated spectra rests on the comparison of the data with 
the predictions of stellar population models. The accuracy 
of such predictions is highly dependent on the input calibra- 
tion of the Ca II triplet line-strengths in terms of the main 
atmospheric stellar parameters, namely effective tempera- 
ture, surface gravity and metallicity. Such calibrations have 
been either theoretical, based on model atmospheres, with 
their associated uncertainties (e.g. Smith & Drake 1987, 
1990; Erdelyi-Mendes & Barbuy 1991; J0rgensen, Carls- 
son & Johnson 1992; Chmielewski 2000), or based on em- 
pirical stellar libraries with a poor coverage of the atmo- 
spheric parameter space (e.g. Jones, Alloin & Jones 1984; 
Diaz, Terlevich & Terlevich 1989; Zhou 1991; Mallik 1994, 
1997; Idiart, Thevenin & de Freitas Pacheco 1997, hereafter 
ITD). For a more detailed discussion of these calibrations 
see Paper I and III. The quality of this calibration has been 
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the major drawback of previous stellar population models 
which have included predictions for the strength of the Ca II 
triplet feature (Vazdekis et al. 1996, hereafter V96; ITD; 
Mayya 1997; Garci'a-Vargas, Molla & Bressan 1998; Molla 
& Garcia- Vargas 2000, hereafter MGV). In this work we 
solve for these deficiencies by updating our model predic- 
tions on the basis of the ample near-IR empirical stellar 
spectral library and calibrations presented in Paper I, II and 
III. Therefore this library constitutes the main ingredient of 
these new models. 

Although this paper is not the first to present mod- 
els for the Can triplet in integrated stellar systems, it is 
the first to present integrated spectra at intermediate reso- 
lution on the basis of an extensive empirical stellar library. 
Bruzual &: Chariot (2003) present spectra in this wavelength 
region based on the theoretical library of Lejeune, Cuisinier 
& Buscr (1997) or the observed library of Pickles (1998), 
both at a resolution of 10-20 A, i.e. an order of magnitude 
larger than our models. On the other hand, Schiavon, Bar- 
buy & Bruzual (2000, hereafter SBB) predicted SSP spectra 
at high resolution employing a theoretical stellar library. 

We must keep in mind that the empirical library used 
here, as well as any empirical library, implicitly includes the 
chemical enrichment history of the solar neighbourhood. For 
instance current models use scaled-solar abundance ratios. 
However there are indications that Ca abundance is not 
enhanced compared to Fe in elliptical galaxies (O'Connell 
1976; Vazdekis et al. 1997, hereafter V97; MGV; Peletier 
et al. 1999; Vazdekis et al. 2001a, hereafter V01A; PS02). 
This discrepancy did mostly show up when investigating the 
Ca 4227 line, which seems to tracks Fe lines in galaxies. Nu- 
cleosynthesis theory (Woosley & Weaver 1995) predicts that 
Ca is an a-element, i.e. is mainly produced in SN Type II 
and therefore should follow Mg. 

The stellar population model that we describe in this 
paper is an improved version, as well as an extension to the 
near-IR spectral range, of the model presented in V96 and 
Vazdekis (1999, hereafter V99). Among the major changes 
introduced here, apart of the stellar library, is the imple- 
mentation of two new IMF-shapes by Kroupa (2001, here- 
after K01) (see § 2.1), the inclusion of the new scaled-solar 
isochrones of the Padova group (Girardi et al. 2000, hereafter 
GOO) (see § 2.2), and the transformation of their theoretical 
parameters to the observational plane (i.e. fluxes and colors) 
on the basis of almost fully empirical photometric stellar li- 
braries, such as the one of Alonso, Arribas & Martinez-Roger 
(1999) (see § 2.3). We also have significantly improved the 
way in which we grid our stellar library for computing a stel- 
lar spectrum of a given set of atmospheric parameters (see 
§ 3). A full description of the behaviour of the Call triplet 
in SSPs of ages larger than 1 Gyr is provided in § 4, whilst 
in § 6 we extend our model predictions to some younger stel- 
lar populations, where the Asymptotic Giant Branch (AGB) 
contribution to the total luminosity in this spectral range is 
very important. In § 5 we discuss the behaviour of the Call 
triplet feature as a function of the spectral resolution and 
galaxy velocity dispersion. Detailed comparisons with previ- 
ous papers are made in § 7. In § 8 we establish the validity of 
these new models by making a direct comparison with galac- 
tic globular clusters and with two galaxies of very different 
luminosities. We use NGC 4478 and NGC 4365, for which 
we have data both in the Can triplet spectral region and 



in the blue (V01A). A detailed analysis of galaxy spectra, 
however, is given in Cenarro et al. (2003, hereafter C03), in 
which we present new Call triplet data for a large sample 
of elliptical galaxies, and show to what extent our view of 
the stellar populations in elliptical galaxies may change as 
a result of the indices in this wavelength region. In § 9 we 
write our conclusions. Finally, we provide three appendices 
with further details on the model construction as well as 
an analysis of the main uncertainties affecting our model 
predictions. 



2 THE MODELS 

The predictions that we present in this paper for the near- 
IR spectral range are calculated on the basis of the evolu- 
tionary stellar population synthesis model presented in V96 
and V99. In this section we describe the main ingredients 
of these models stressing those updates that have been per- 
formed since the original version. In § 2.1 we summarize 
the IMF types used by our models and introduce two new 
IMF shapes considered in this paper. In § 2.2 we describe 
the properties of the new isochrones of GOO that are em- 
ployed here, showing the relevant differences with respect to 
the ones used in V96 (mostly Bertelli et al. 1994, hereafter 
B94). In § 2.3 we describe how we transform the theoretical 
parameters of these isochrones to the observational plane, 
which is performed on the basis of extensive empirical stel- 
lar libraries rather than using model stellar atmospheres. 
Finally, in § 2.4 we provide a brief summary of the main 
predictions of our models. 

2.1 The initial mass function 

In this paper we adopt the IMF shapes described in V96 
(i.e unimodal and bimodal) and the two IMFs recently in- 
troduced by K01. To calculate the number of stars in the 
mass interval mtom + dm we use the following analytical 
approaches for the IMF: 

• Unimodal: a power-law function characterized by its 
slope /i as a free parameter 

$(m) ocm- (M+1) . (1) 

It is worth noting that the Salpeter (1955) IMF corresponds 
to /! = 1.3. 

• Bimodal: similar to the unimodal IMF for stars with 
masses above 0.6 Mq, but decreasing the number of the stars 
with lower masses by means of a transition to a shallower 
slope (see V96) . This IMF is also characterized by the slope 
ti. 

• Kroupa (2001) universal: a multi-part power- law IMF, 
which is similar to the Salpeter (1955) IMF for stars masses 
above 0.5 M©, but with a decreasing contribution of lower 
masses by means of two flatter segments. 

• Kroupa (2001) revised: a multi-part power-law IMF, in 
which the systematic effects due to unresolved binaries on 
the single-star IMF have been taken into account. This IMF 
is steeper than the universal IMF by A/i ~ 0.5 in the mass 
range 0.08 Mq < m < 1 Mq. 

These IMF types are plotted in Figure 1. Further details 
of the IMF definitions are given in Appendix A. We also 
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Figure 1. IMF types used by our models. Although not 
here the slope of the power-law (p) is allowed to vary 
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refer the reader to V96 and KOI for a full description of 
these IMFs. 



2.2 Isochrones 

The new model makes use of the updated version of the theo- 
retical isochrones of the Padova group (GOO), whereas in V96 
and V99 we used the B94 isochrones and the stellar tracks 
of Pols et al. (1995) for very low-mass stars (M < O.6M ). 
These isochrones cover a wide range of ages and metallic- 
ities and include the latest stages of the stellar evolution 
through the thermally pulsing AGB regime to the point 
of complete envelope ejection (employing a synthetic pre- 
scription). The GOO set has a larger age resolution for ages 
lower than 10 Gyr, however the largest metallicity covered 
is Z = 0.03 (instead of Z = 0.05 as in B94). Hereafter we 
will refer to the metallicity following the relation 



[M/H] = log(Z/Z ), 



(2) 



where Zq — 0.019 according to the reference value adopted 
by GOO. The lower mass cutoff of the new isochrones has 
been extended down to O.15M . The input physics of the 
isochrones have been updated with an improved version of 
the equation of state, the opacities of Alexander & Ferguson 
(1994) (which make the Red Giant Branch, RGB, slightly 
hotter than in B94 isochrones, see Figure C2) and a milder 
convective overshoot scheme. A helium fraction was adopted 
according to the relation: Y ~ 0.23 + 2.25Z. 

It is worth noting that, although giant elliptical galax- 
ies show non-solar abundance trends (e.g. Peletier 1989; 
Worthey, Faber & Gonzalez 1992; K00; V01A), only scaled- 
solar abundance ratios are adopted for this set of isochrones. 
In particular, these galaxies seem to show an enhancement 



of the a-elements (e.g. O, Ne, Mg, Si, S, etc.) with respect 
to Fe. For a given total metallicity, the a-enhanced mix- 
tures yield lower opacities, which translates into an increase 
of the temperature of the stars in the Main Sequence (MS) 
and the RGB phases (e.g. Salaris & Weiss 1998, hereafter 
SW98; VandenBerg et al. 2000; Salasnich et al. 2000, here- 
after S00; Kim et al. 2002). Atomic diffusion is not included 
in the present set of isochrones, although this phenomenon 
decreases the temperatures of the stars of the turnoff (e.g. 
Salaris, Groenewegen & Weiss 2000), an effect which has 
been found to be useful to be able to fit the low line-strength 
values of the Balmer lines in the spectra of metal-rich glob- 
ular clusters (Vazdekis et al. 2001b, hereafter V01B). 



2.3 Transformation to the observational plane 

To transform the theoretical parameters of the isochrones 
to the observational plane, i.e. colours and magnitudes, we 
make use of relations inferred on the basis of extensive em- 
pirical photometric stellar libraries (rather than by imple- 
menting theoretical stellar atmospheric spectra) to obtain 
each colour as a function of the temperature, metallicity 
and gravity. We use the metallicity-dependent empirical re- 
lations of Alonso, Arribas & Martmez-Roger (1996, 1999), 
for dwarfs and giants respectively (each stellar sample is 
composed of ~ 500 stars) . The derived temperature scale is 
based on the IR-Flux method and, therefore, is marginally 
dependent on the model atmospheres. This treatment for the 
giants is the most important difference with respect to the 
models of V96, where we used the empirical calibration of 
Ridgway et al. (1980) to obtain the V — K colour from T e g , 
after which we applied the colour-colour conversions of John- 
son (1966). The empirical (not the theoretical) compilation 
of Lejeune, Cuisinier & Buser (1997, 1998) (and references 
therein) are used for the coolest dwarfs (T c g < 4000 K) and 
giants (T ff < 3500 K), respectively, for solar metallicity; a 
semi-empirical approach was applied to other metallicities 
on the basis of these relations and the model atmospheres 
of Bessell et al. (1989, 1991) and the library of Fluks et al. 
(1994). The empirical compilation of Lejeune et al. was also 
used for stars with temperatures above ~ 8000 K. 

Finally, we use the metal-dependent bolometric cor- 
rection given by Alonso, Arribas & Martmez-Roger (1995, 
1999) for dwarfs and giants, respectively. For the Sun we 
adopt the bolometric correction BCq = —0.12, with a bolo- 
metric magnitude of 4.70 according to these authors. 



2.4 Summary of the previous model predictions 

Using data obtained at different spectral ranges, or com- 
bining spectroscopic and photometric information, increases 
our constraining power for interpreting the stellar popula- 
tions (e.g. V97). Furthermore, the use of a set of model pre- 
dictions for a given spectral range, such as the ones pre- 
sented in this paper and those of other authors for a differ- 
ent spectral range, might drive to systematic effects due to 
differences in the adopted model prescriptions (e.g. Chariot, 
Worthey & Bressan 1996). Furthermore, it is well known 
that the ages or metallicities inferred by means of stellar 
population synthesis modeling should be taken on a relative 
basis (e.g. V01B). Therefore we briefly summarize below our 



© 2000 RAS, MNRAS 000, 1-27 



The nIR Call triplet: stellar population synthesis models 5 



model predictions for those interested in applying the mod- 
els presented in this paper. 

Broadband colours, mass-to-light ratios and line- 
strengths for the main features in the Lick/IDS system 
were predicted in V96. A new version of these models in- 
clude the H<5 and H7 indices of W097, as well as the 
break at -4000 A of Gorgas et al. (1999). SSP spec- 
tral energy distributions in the blue and visual at res- 
olution 1.8 A (FWHM) were presented in V99. New H7 
index predictions can be found in VA99 and V01B. Fi- 
nally, optical and near-IR surface brightness fluctuations 
magnitudes and colours (including the WFPC2-HST fil- 
ter system) were presented in Blakeslee, Vazdekis & Ajhar 
(2001), where we include a description of the updated mod- 
els. All these predictions are available at the web page: 
http : //www . iac . es/galeria/ vazdekis/. 



3 NEW MODEL PREDICTIONS FOR THE 
NEAR-INFRARED 

Using the model ingredients described in § 2 and the new 
stellar spectral library presented in Paper I and II, we fol- 
low two different approaches two build up two set of model 
predictions. First we compute spectral energy distributions 
at resolution (FWHM=1.5A) for different metallicities, ages 
and IMFs. In § 3.1 we further the preparation of our stellar 
library for stellar population synthesis modeling, whereas 
in § 3.2 and Appendix B we describe our model and com- 
puting details. In Appendix C we describe some tests that 
we have performed to show the main uncertainties affect- 
ing the obtained model predictions, as well as to probe the 
robustness of our approach. For the second set of model 
predictions, which is described in § 3.3, we make use of the 
empirical fitting functions presented in Paper III (based on 
the same stellar spectral library) to calculate the strengths 
of the Ca 11 triplet feature by means of the index definitions 
given in Paper I. Finally, in § 3.4 we summarize these two 
sets of model predictions and describe the spectral proper- 
ties and parameter coverage of the synthesized SSP spectral 
energy distributions. 



3.1 The near-IR empirical stellar spectral library 

The main ingredient of the models presented in this paper is 
a new stellar library covering the spectral range around the 
near-IR Call triplet feature at resolution 1.5 A (FWHM). 
A full description of this library was given in Paper I, where 
we also propose a new index definition for this feature (i.e. 
CaT*). In order to make this stellar library useful for stellar 
population synthesis modeling, we require an homogeneous 
set of accurate atmospheric parameters (Toff, log g, [Fe/H]) 
for the stellar sample (see Paper II). As pointed out in V99, 
this step is particularly important to avoid systematic trends 
in the parameters among different authors. In this section 
we give further details on the preparation of the stellar li- 
brary for its implementation in the model. This requires us 
to identify those stars whose spectra might not be properly 
representing a given set of atmospheric parameters (§ 3.1.1). 
Another important step is the characterization of the pa- 
rameter coverage of the stellar library. This step, which is 



described in § 3.1.2, allows us to to understand the limita- 
tions of the models presented in § 3. 

3.1.1 Preparation of the stellar library 

To optimize the stellar library we have performed a second 
selection on the stars to use for the stellar population mod- 
eling. For this purpose we identified all the spectroscopic 
binaries, making use of the SIMBAD database, as well as 
those stars with high signal of variability (AV > 0.10 mag). 
We used the Combined General Catalogue of Variable Stars 
of Kholopov et al. (1998) (the electronically-readable ver- 
sion provided at CDS). Moreover, we removed from the 
sample a number of these stars that were found irrelevant 
since the stellar library contains a reasonable large number 
of stars with similar atmospheric parameters with no such 
high variability feature. The removed stars were: HD 10476, 
HD 36079, HD 113139, HD 113226, HD 125454, HD 138481, 
HD 153210 and HD 205435. We refer the reader to § 3.2.2 
and Appendix B for a full description of how we deal with 
these stars when calculating a representative stellar spec- 
trum for a given set of atmospheric parameters. 

A number of stars were removed due to the high resid- 
uals obtained when calculating the Ca n triplet fitting func- 
tions (see Paper III). These stars were: HD 17491, HD 35601, 
HD 42475, HD 115604, HD 120933, HD 121447, HD 138279, 
HD 181615, HD 217476 and HD 222107. 

We did remove from the original stellar sample most 
of the stars for which at least one of the three atmospheric 
parameters was lacking (see Table 5 of Paper II) . However we 
kept 47 of these stars, with unknown metallicities, but whose 
temperatures were either larger than 9000 K or smaller than 
4000 K. 

We are interested in achieving the largest possi- 
ble spectral resolution for our model predictions, i.e. 
FWHM = 1.5 A. We therefore removed from the library the 
few stars observed with instrumental configurations provid- 
ing — 2.1 A, as described in Table 2 of Paper I (i.e. those 
corresponding to the observing runs 5 and 6). These stars 
were those corresponding to NGC 188 as well as NGC 7789 
676, NGC 7789 859, NGC 7789 875, NGC 7789 897 and 
NGC 7789 971. The remaining stars whose spectra were ob- 
tained at resolutions slightly better than 1.5 A were broad- 
ened to match this resolution (see Paper I). 

For different technical reasons we removed from the 
sample an additional number of stars. In particular, a sub- 
set of stars were discarded because their spectra were of low 
signal-to-noise (i.e. M 5 11-53, M 92 1-10 and M 92 XII-24). 
M 71 1-63 and HD 232979 were removed for being affected 
by cosmic rays and M 67 F119 for the poor flux calibration 
quality. 

As a result of this cleaning process we obtained a sub- 
sample composed of 616 stars. A final step was to revise the 
atmospheric parameters presented in Paper II for a num- 
ber of stars. This is motivated by the evidence that their 
line-strength values significantly deviate from the expected 
general trends predicted by the fitting functions. In partic- 
ular, we modified the T cff of HD 167006 from 3470 K to 
3650 K in order to bring it in better agreement with the 
value predicted by a new set of fitting functions, which cal- 
culates the slope of the spectrum around the Can triplet 
feature due to the TiO molecular bands (Cenarro 2002). Fi- 
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Figure 2. The fundamental parameter coverage of the subsamplc 
of selected stars 



nally, we also changed the metallicity adopted for the stars 
of M 71 from -0.70 (Carretta & Gratton 1997) to -0.84 
in order to achieve a better agreement with the predictions 
of the fitting functions presented in Paper III for the Call 
triplet feature. As described in § 3.3 of Paper III, the stars 
of this globular cluster showed systematic residuals in com- 
parison to the strengths predicted by the fitting functions 
(i.e. ACaT* = —0.35), which translates to the metallicity 
value that we adopt here. We refer the reader to that paper 
for an extensive discussion on the subject. 



3.1.2 Stellar atmospheric parameter coverage 

Figure 2 shows the parameter coverage of the selected sub- 
sample of stars. The adopted values for the T e tt, logg and 
[Fe/H] of the stars of our sample were extensively discussed 
in Paper II. This figure includes the corrections described in 
§ 3.1.1. The figure shows that all types of stars are well rep- 
resented for solar metallicity. We also see a good coverage 
of typical dwarfs and giants of 5000 < T e g < 6500 K and 
4000 < T c ff < 5500 K respectively, for metallicities in the 
range -0.7 < [Fe/H] < +0.2. 

Unfortunately, very metal rich stars, particularly 
dwarfs, are scarce in the sample preventing us to provide 
predictions for SSPs of metallicities larger than +0.2. We 
also note that very cool dwarfs and giants of non solar 
metallicity are scarce. These dwarfs are particularly use- 
ful for predicting stellar populations with very steep IMFs, 
i.e. dwarf-dominated, whereas the effect of these low metal- 
licity giants will be discussed below. An important gap is 
the absence of metal-poor dwarfs of T e s > 6500 K. The 
lack of these stars do not allow us to predict SSPs of ages 
younger than ~6 Gyr for [Fe/H] ~ —0.7 and younger than 
~ 2.5 Gyr for [Fe/H] ~ —0.4, since they will be representing 



the predicted Turnoff for these stellar populations. Over- 
all we conclude that the metallicity range where we are 
mostly safe is —0.7 < [Fe/H] < +0.2. With more caveats 
we also can predict stellar populations for lower metallici- 
ties where the contribution of dwarf stars with temperatures 
above T e g > 6500 K is not significant, and the Horizontal 
Branch is not too blue (note in Fig. 2 the lack of metal-poor 
giants with temperatures larger than ~ 6000 K). Therefore 
in this metallicity range we mostly need to focus on SSPs 
with ages larger than ~ 10 Gyr, but smaller than ~ 13 Gyr. 



3.2 Spectral synthesis of SSPs 

Predicting spectral energy distributions for stellar popula- 
tions at intermediate or moderately high resolution, rather 
than the strengths of a given number of features (see § 3.3), 
requires a complete stellar spectral library where all type 
of stars are well represented, and whose spectra are of high 
quality and flux calibrated. Theoretical libraries are difficult 
to achieve due to the limitations of the input physics and the 
computing power. On the other hand, the empirical libraries 
actually suffer from the difficulties in observing a good num- 
ber of representative stars of all types whose atmospheric 
parameters are accurately known. Spectral energy distribu- 
tions at very low resolutions for stellar populations covering 
a large range of ages and metallicities were predicted, ei- 
ther based on the Kurucz (1992) theoretical stellar library 
(e.g. Bressan, Chiosi & Fagotto 1994; Kodama & Arimoto 
1997; Kurth, Fritze-V.Alvensleben & Fricke 1999; Poggianti 
ct al. 2001) or empirical stellar libraries (Bruzual & Chariot 
1993). However the low dispersion of the predicted spectra 
did not allow to measure reliable line-strengths. 

More recently, the availability of such libraries at higher 
dispersion has made it possible to follow the spectral synthe- 
sis approach for the blue spectral range. In V99, we used the 
extensive empirical stellar library of Jones (1999), composed 
of more than 600 stars, to predict spectra of SSPs at reso- 
lution 1.8A (FWHM). This approach allowed us to analyse 
the entire spectrum of a galaxy at one time and to measure 
line-strengths at the resolution imposed by the instrumental 
configuration of the observational data, or to galaxy inter- 
nal velocity dispersion. This represents a great advantage 
over previous models, such as those of Worthey (1994) or 
V96, which only predicted the strengths of a given number 
of features under a determined instrumental configuration 
and resolution, i.e. the Lick/IDS system (Gorgas et al. 1993; 
W94). Note that the old models required the observers to 
transform their data to the characteristics of the model pre- 
dictions. A different approach was followed by SBB, who 
synthesized SEDs of SSPs for the near-IR spectral range on 
the basis of a fully synthetic stellar spectral library by SBB. 
In the following, we describe in detail how we make use of 
our stellar library for predicting near-IR spectra at resolu- 
tion 1.5A (FWHM) for SSPs of different ages, metallicities 
and IMF slopes. Although we have followed a procedure 
similar to that of V99 for the optical spectral range we have 
updated the method as described below. 

To predict any integrated observable (e.g. colours, line- 
strengths, SEDs) of an SSP of age t and metallicity [M/H] 
the synthesis code yields a stellar distribution on the basis of 
the adopted isochrones and IMF (see V96 for more details) . 
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Therefore, the integrated spectrum of the SSP, S\(t, [M/H]), 
is calculated in the following way: 

S x (t,[M/H}) = / Sx(m,t,[M/U])N(m,t) x 

J m\ 

F AAref (m,t,[M/H])dm, (3) 

where S\(m, t, [M/H]) is the empirical spectrum correspond- 
ing to a star of mass m and metallicity [M/H] which 
is alive at the age assumed for the stellar population t. 
^AA ref (jn,t, [M/H]) is its corresponding absolute flux at a 
certain wavelength reference interval, AA ro f, and N(m,t) 
is the number of this type of stars, which depends on the 
adopted IMF. mi and m t are the stars with the smallest 
and largest stellar masses, respectively, which are alive in 
the SSP. The upper mass limit depends on the age of the 
stellar population. 

The spectrum to be assigned to each of the stars 
required by Eq. 3 is selected from the empirical stellar 
database, whereas the corresponding absolute flux is as- 
signed following the prescriptions of our code, which are 
based on the extensive empirical photometric stellar libraries 
described in § 2.3. In § 3.2.1 and § 3.2.2 we describe how we 
calculate i<AA ref (m, t, [M/H]) and S\(m, t, [M/H]), respec- 
tively. 

3.2.1 Assigned stellar flux 

We need to find a common reference wavelength interval, 
AA rc f, from which we can scale S\{m, t, [M/H]). We selected 
for this purpose the spectral range AA8475-8807A around 
the peak of the I filter of Johnson (1966). We have chosen 
the Johnson (1966) I filter because its effective wavelength 
is closer to the spectral range covered by our stellar spectra 
than the I filter of the Cousins system (Bessell 1979). We 
divide our stellar spectra by the average flux in the selected 
wavelength reference interval. To calculate the absolute flux 
in this spectral region for each of the stars requested by the 
code, -FXA ref (fn, t, [M/H]), we first used equation (6) of Code 
ct al. (1976) which provides a relation between the absolute 
V flux and the calibrated bolometric correction in order to 
infer equation (2) of V99 from which we derive the following 
relation: 

F, ln +o.4(v-/) 
1 = — _ U) 

F holQ 10+0-4(M v -3.762) ' V > 

where Fi is the absolute flux in the Johnson I-band. The 
next step is to calculate the absolute flux per angstrom in 
the selected reference interval: 

where /„ represents the fraction of Fi corresponding to the 
normalization wavelength interval, i.e. A A re f =8475-8807 A. 
This factor is calculated on the basis of the empirical stellar 
spectral library of Pickles (1998), for which we have com- 
pared the flux in the I filter and in the normalization region, 
obtaining the following relations 

/„=-1.620675 + 0.496916 log T eff , log T eS < 3.55,dwarfs 
/„=-0.163463 + 0.086852 log T cS , log T eS < 3.55,giants (6) 
/„=+0.156952 - 0.003324 log Teff, log T eff > 3.55 . 



0.1 60 F 



0.1 50 F 




0.100 F I : 

3.4 3.6 3.8 4.0 4.2 4.4 4.6 
l °9 T eff 

Figure 3. Normalization wavelength interval flux correcting fac- 
tor versus logT e ff. The fit is given in Eq. 6. 



Note that f n is nearly constant, except for M stars with 
temperatures lower than ~ 3500 K. See Figure 3. 

Figure 4 illustrates the contributions of the different 
evolutionary phases to S\(t, [M/H]) in the selected normal- 
ization wavelength interval, according to the values calcu- 
lated for F&\ re , for each star. The figure shows the fraction 
of luminosity (in percentage) of these evolutionary stages 
for SSPs of several IMF slopes (unimodal), metallicities and 
ages. The RGB, MS and HB are the main contributors to 
this spectral range. The weight of the AGB decays very 
rapidly as the stellar population evolves, whereas the RGB 
is completely build-up when the age reaches ~ 2 Gyr. The 
contribution of the Sub Giant Branch (SGB) is the small- 
est one in all the diagrams. The effect of the IMF is mostly 
reflected in the relative contribution of the MS and RGB. 

Finally, it is worth recalling that the fact that the ef- 
fective wavelength of the Johnson I filter experiments a sig- 
nificant shift as a function of the spectral type, makes our 
scaling approach more secure than if we had chosen, as a 
reference interval, a significantly narrower spectral range. 
In fact we have tested to scale the stellar spectra with Fj 
(i.e. non corrected by /„), selecting as a normalization region 
either the pseudocontinuum AA8474-8484A (i.e. cl, see Ta- 
ble 1) or 8776-8792A (i.e. c5). For the synthesized SSP spec- 
trum, e.g. for 12.6 Gyr and solar metallicity, the slope of the 
continuum in the wavelength range AA8500-8850A (mea- 
sured as Flux C 5/Flux c i) varies by ~ 4%, and the value of 
the CaT* index by ~ 0.14A (which represents ~2% of the 
index strength) . This effect is due to the fact that the selec- 
tion of a redder continuum emphasizes the contribution of 
redder stars in Eq. 3. 
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Figure 4. Fraction of luminosity (in percentage) of the different stellar evolutionary phases in the Ca II triplet spectral region for SSPs 
of unimodal IMF of different slopes (from top to bottom panels), metallicities (from left to right panels) and ages (in each panel). 
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3.2.2 Adopted stellar spectrum 

In Appendix B we describe the method that we have followed 
for obtaining from the whole subsample of selected stars, a 
representative spectrum, S\{m,t, [M/H]), for a given set of 
atmospheric parameters 9, \ogg and [M/H]. It is worth not- 
ing that rather than T c ff we have chosen to work with 9 
(= 5040/T e ff) in order to follow the scale adopted in Pa- 
per III for calculating the empirical fitting functions for the 
Ca II triplet feature. We have tested however that using this 
parameter does not affect the resulting integrated SSP spec- 
trum. 

Basically our approach consists in finding all the stars 
whose parameters are enclosed within a given box in the 
3-parametrical space. The method ensures the presence of 
stars in several directions to alleviate the effect of asym- 
metries in the distribution of stars around the point. For 
example, for metallicity larger than solar, more stars are 
usually found with lower, rather than higher, metallicity. 
The size of the selected finding box is taken to be a function 
of the density of stars around the requested point, i.e. the 
larger the density the smaller the box. Moreover, in the more 
populated parametrical regions, the typical uncertainties in 
the determination of the atmospheric parameters are usually 
smaller than in regions with lower densities (see Paper II), 
in agreement with the adopted criterion. When needed, the 
starting box is enlarged until finding suitable stars. 

The obtained stars are combined according to their 
parameters 9, log<? and [M/H] and to the signal-to-noise 
of their spectra. However, the method ensures that spec- 
troscopic binaries and stars with anomalous signatures of 
variability do not significantly contribute to the average 
spectrum. Finally, we apply a minimal correction to these 
weights in order to ensure the obtention of a spectrum whose 
atmospheric parameters matches the requested ones. 



3.3 Line-strength index synthesis 

The most widely employed approach for modeling and 
studying stellar populations using absorption lines is based 
on the, so-called, empirical fitting functions. These functions 
describe the strength of, previously defined, spectral features 
in terms of the main atmospheric parameters. These calibra- 
tions are directly implemented into the stellar population 
models to derive the index values for stellar populations of 
different ages, metallicities and IMFs. The main advantage 
of this approach over the spectral synthesis, which provides 
full spectra, is that the later method relies on the availability 
of a complete stellar library at the appropriate spectral reso- 
lution. The difficulties in achieving such libraries in the most 
recent past made the fitting functions method the only pos- 
sible approach to provide robust diagnostics for the stellar 
populations (e.g. Buzzoni 1993; Worthey 1994; V96; Tan- 
talo, Chiosi & Bressan 1998; Kurth et al. 1999; Maraston 
& Thomas 2000; Poggianti et al. 2001; Bruzual & Char- 
lot 2003). The predictions based on these models have been 
widely used along the last decade (see Trager et al. 1998 for 
an extensive review). 

In the blue and visible part of the spectrum, we pre- 
dicted in V96 the most important absorption line-strengths 
of the Lick/IDS system, at intermediate resolution FWHM 
~ 9 A. We used the empirical fitting functions of Gorgas et 



Table 1. Bandpass limits for the generic indices CaT and PaT. 



CaT central PaT central Continuum 

bandpasscs (A) bandpasses (A) bandpasses (A) 

Cal 8484.0-8513.0 Pal 8461.0-8474.0 cl 8474.0-8484.0 
Ca2 8522.0-8562.0 Pa2 8577.0-8619.0 c2 8563.0-8577.0 
Ca3 8642.0-8682.0 Pa3 8730.0-8772.0 c3 8619.0-8642.0 

c4 8700.0-8725.0 
c5 8776.0-8792.0 



al. (1993), W94 and W097. Moreover we also predicted in 
that paper the strengths of the Call triplet on the basis of 
the stellar spectral library and the index definition of Diaz 
et al. (1989). Call triplet strengths were also calculated by 
Garcia- Vargas et al. (1998), MGV and ITD (see § 7 for a 
full comparison with these models). It is worth noting that 
the use of these predictions requires one to adapt the data 
to the characteristics of the instrumental-dependent stellar 
spectral library employed by these models (see W097 for an 
extensive review of the method). 

Since we now have a better stellar library and models, 
our previous predictions for the Ca II triplet presented in V96 
are superseeded by the ones calculated in this section. Here 
we adopt for the Can triplet feature the index definition 
given in Paper I (CaT*). This new index removes the H 
Paschen lines contamination and it is, therefore, a reliable 
indicator of pure Call triplet strength. The index is given by 
CaT* = CaT - 0.93 PaT, where the CaT and PaT indices 
measure the strengths of the raw calcium triplet and of three 
pure H Paschen lines, respectively (see Table 1). All these 
indices are computed for a nominal resolution FWHM = 
1.5A. Further details, as well as relations to transform this 
index to other popular definitions for the Ca 11 triplet feature 
can be found in Paper I. 

Empirical fitting functions for the CaT and PaT indices 
were calculated in Paper III. These functions are polynomia 
that relate the strength of these indices to the stellar atmo- 
spheric parameters (9, log g, [M/H]) derived in Paper II. The 
new fitting functions reveal a complex behaviour of the Ca 11 
triplet features as a function of the atmospheric parameters. 
In particular, for hot and cold stars, the temperature and 
luminosity class are the main driving parameters, whereas, 
in the mid-temperature regime, the three atmospheric pa- 
rameters play an important role. Among the advantages of 
the new fitting functions over previous predictions is that 
they include the whole range of effective temperatures and, 
in particular, cold stars. We refer the interested reader to 
Paper III for a full description of these functions. 

The implementation of these fitting functions in our 
stellar population code is readily done as follows. A flux- 
weighted index i for a SSP of age (t) and metallicity ([M/H]) 
is given by: 

J*™* i(m,t, [M/B])N(m,t)F A x Io{ (m,t, [M/H])dm 
tSSP ^ Q N(m,t)F AKaf (m,t,[M/U])dm ' (?) 

where the flux i*AA ref (m, t, [M/H]) is used to scale the contri- 
bution of a star of mass m and age t in a given evolutionary 
stage. We adopt the flux corresponding to the scaling wave- 
length interval used for the spectral synthesis, i.e. AA re f, 
since it is centered on these features. 
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3.4 Summary of the new model predictions for 
the near-IR 

We summarize in this section the new model predictions 
presented here for the near-IR spectral range. There are two 
main products interesting for the stellar population anal- 
yses: i) the spectral library for SSPs, for which we list in 
Table 2 the main properties of the synthesized SSP spectra, 
as well as the model parameters coverage according to the 
limitations described in § 3.1.2, and ii) the CaT*, CaT and 
PaT index strengths predicted on the basis of the empiri- 
cal fitting functions, which are listed in Table 3. It is worth 
noting that the index values given in this table for SSPs of 
different ages, metallicities and IMF types and slopes, arc 
calculated at resolution 1.5 A (FWHM). Therefore, users 
interested in applying these predictions should correct their 
index measurements to the ones they would have obtained 
at the resolution of the models presented here (see the rela- 
tions given in § 5). 

These two set of model predictions can either be ob- 
tained from the authors or from the following web pages: 
http : //www .ucm . es/inf o/Astrof / ellipt/CATRIPLET . html, 
http : //www . nott ingham . ac . uk/~ppzrf p/CATRIPLET . html , 
and http : / /www . iac . es/ galeria/ vazdekis/. 



4 THE BEHAVIOR OF THE CAII TRIPLET 
FEATURE IN SSPS 

In this section we focus on describing the behaviour of the 
integrated Call triplet feature as a function of the relevant 
stellar population parameters in old-aged SSPs. In § 6 wc 
will present a discussion for some younger ages. We also dis- 
cuss here some of the main differences seen in the full spectra 
predicted by the models, but we refer the reader to C03 and 
Cenarro (2002), for a detailed description of the behaviour of 
other features and the slope of the characteristic continuum 
around the Call triplet feature, which is largely determined 
by the contributions of the TiO molecular bands. In § 4.1 we 
describe the time evolution of the Call triplet feature and 
its dependence on metallicity, whereas in § 4.2 we discuss 
the effects of the IMF. 

4.1 The effects of the age and metallicity 

In Figure 5 we plot the time evolution of the CaT* index 
as calculated on the basis of the fittings functions of Pa- 
per III (see § 3.3) and as measured on the synthesized spec- 
tra (§ 3.2). The two set of predictions are given at resolution 
1.5A (FWHM), i.e. a = 22.2 km s _1 (the behaviour of the 
Call triplet feature as a function of spectral resolution will 
be discussed in § 5). Fig. 5 reveals three major features: i) 
the strength of the CaT* index does not vary much for stel- 
lar populations older than ~ 3 Gyr, for all metallicities, ii) 
there is a strong sensitivity of the CaT* index as a func- 
tion of metallicity for values below [M/H] ~ —0.5 and Hi) 
for higher metallicities, the CaT* index does not show any 
significant dependence on the metallicity and the strength 
of this feature tends to saturate. 

The observed saturation of the CaT* index can be un- 
derstood in terms of comparing the contributions of the MS 
and RGB to the total luminosity of an SSP in this spectral 
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Figure 5. The CaT* index in SSPs as a function of age and 
metallicity. Different line types refer to the metallicity of the SSP, 
as quoted in the plot. All the models are calculated with a uni- 
modal IMF of slope fi = 1.3. We plot the CaT* index as calculated 
on the basis of the fittings functions of Paper III (ff, thin lines), 
and as measured on the synthesized spectra (ss, thick lines). The 
spectral resolution is 1.5A (FWHM) (i.e. a = 22.2 km s" 1 ) for 
the two set of predictions. 



range (see Fig. 4). If we focus on the second row of panels 
of this figure, i.e. for [i = 1.3 (we will discuss the effects 
of the IMF in § 4.2), we see that the relative contribution 
of the MS with respect to the RGB gets larger as metallic- 
ity increases. Dwarf stars provide smaller strengths for this 
feature than giants according to the results obtained in Pa- 
per III. Moreover, in metal-rich SSPs the RGB is populated 
with an increased fraction of very cool stars. For tempera- 
tures cooler than ~ 3500 K, the CaT* index does not de- 
pend on metallicity and decreases very rapidly with decreas- 
ing temperature (see Fig. 7 of Paper III). The saturation of 
the Ca II triplet for large metallicities is a result that differs 
from the predictions of previous authors (see § 7). This re- 
sult might be supported by the CaT* index measurements 
obtained for a large sample of ellipticals (C03) and bulges 
of spirals (Falcon-Barroso et al. 2002), as well as for the 
most metal- rich galactic globular clusters (see § 8.1). Fur- 
thermore, Cohen (1979), Bica & Alloin (1987), Terlevich et 
al. (1990b) and Houdashelt (1995) found that the strength of 
the Ca II triplet did not vary much among early-type galax- 
ies of different absolute magnitudes and colours. 

Fig. 5 shows a reasonably good agreement between our 
CaT* index predictions based on the empirical fitting func- 
tions and the ones measured on the synthesized spectra. 
We note, however, a better agreement for larger metallici- 
ties, where the CaT* strengths provided by the spectral syn- 
thesis are marginally larger than the ones predicted on the 
basis of the empirical fitting functions. For [M/H] = —0.38 
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Table 2. Spectral properties and parameter coverage of the synthesized model SEDs 



Spectral properties 


Spectral range 


AA8348.85-8950.65 A 


Spectral resolution 


FWHM = 1.5A, cr = 22.2 km s" 1 


Linear dispersion 


0.85 A/pix (29.358 km s" 1 ) 


Continuum shape 


Flux-scaled 


SSPs parameter coverage 


IMF type 


Unimodal, Bimodal, Kroupa universal, Kroupa revised 


IMF slope (for unimodal and bimodal) 0.3 - 3.3 


Mctallicity 


-1.68, -1.28, -0.68, -0.38, 0.0, +0.20 


Age ([M/H] < -0.68) 


10.0 < t < 17.78 Gyr 


Age ([M/H] = -0.68) 


5.62 < t < 17.78 Gyr 


Age ([M/H] = -0.38) 


2.51 < t < 17.78 Gyr 


Age ([M/H] = 0.0) 


0.1 < t < 17.78 Gyr 


Age ([M/H] = +0.2) 


1.0 < t < 17.78 Gyr 



and [M/H]= —0.68, the spectral synthesis approach provides 
smaller CaT* values, and the obtained differences between 
the predictions based on the two approaches can be as large 
as 0.2 A (i.e. ~ 3% of the CaT* index value) for [M/H] 
= -0.68. Also, the SSP spectra for [M/H] = -1.28 and 
— 1.68 provide smaller CaT* strengths than the values ob- 
tained on the basis of the fitting functions, although the ob- 
served differences are smaller than those obtained for [M/H] 
= —0.68. These differences must be attributed to differences 
in the methods employed by the two approaches to provide, 
for a given set of stellar atmospheric parameters, a repre- 
sentative average spectrum and index value, respectively. 
One of the main differences between these two methods is 
that the spectral synthesis works at much higher resolution 
in the stellar parameter space than the one used for calcu- 
lating the fitting functions. In the boxes-method employed 
in the spectral synthesis those stars whose parameters re- 
semble more the requested ones are assigned significantly 
larger weights when computing a representative spectrum 
(see Appendix B). Therefore, since the spectral synthesis 
emphasizes the contribution of the most likely stars, the ob- 
tained differences between the two set of model predictions 
are expected to be larger when a given set of stars, with 
a relevant contribution within a stellar population, is lack- 
ing. In fact the lack of giant stars of temperatures lower 
than ~ 4200 K (see Fig. 2) causes the main differences 
seen for [M/H] = -0.38 and -0.68. However, for [M/H] 
< —1.0 these differences are smaller because the coolest stars 
along the isochrones have larger temperatures than those for 
more metal-rich stellar populations and, therefore, the lack 
of these stars constitutes a less severe problem. However, 
we note that these differences in the predictions of the two 
model approaches are of the order of typical observational er- 
rors in galaxy spectra (see Paper I and C03). Moreover, these 
differences are considerably smaller than the ones obtained 
from comparisons between different author predictions (see 
§7). 

It is worth noting that the SSP spectra obtained for 
[M/H] = -1.28 and [M/H] = -1.68 are less reliable than 
the ones for higher metallicities due to the lack, in our stel- 
lar library, of spectra representing important evolutionary 
phases. In particular, this is the case for stars correspond- 
ing to the MS turnoff (i.e. dwarfs with temperatures larger 



than 6000 K), or the Blue Horizontal Branch, which is more 
prominent for SSPs of ages above 13 Gyr (see Fig. 2). This 
is the reason why we only show in Fig. 5 the CaT* values 
measured in the spectra corresponding to SSPs for ages in 
the range 10-13 Gyr, where the contribution of these stars 
is smaller. 

Figure 6 shows the effects of the metallicity (top pan- 
els), age (middle) and IMF (bottom) on the synthesized 
SSP spectra at resolution a = 200 km s _1 (left panels) and 
a — 22.2 km s _1 (right panels). We overplot in the upper 
panels various representative SSP spectra of similar age (i.e. 
10 Gyr) and IMF (unimodal, fi = 1.3), and with different 
metallicities. We have chosen to normalize all the spectra to 
the continuum at AA8619-8642A (see Table 1), to be able 
to see the effect of the metallicity on the third line of the 
Can triplet feature, which can be taken as representative 
of the overall Can triplet feature. We see that the depth of 
this line is virtually constant in the metallicity range (i.e. 
-0.7 < [M/H] < 0.2). We obtain similar results for the 
other lines if normalizing to the appropriate continua. We 
note, however, a strong variation of the slope of the contin- 
uum around the Can triplet feature. For a given age and 
IMF, the larger the metallicity the larger the slope of this 
continuum. The overall shape of the continuum for SSPs of 
high metallicities shows strong similarities to those of early 
M-type stars. In this spectral range the contribution of these 
cool stars to the total luminosity is significantly larger (see 
Fig. 4) than in the visible, where the integrated spectra re- 
semble those of K-giants. The slope of these spectra and its 
comparison to galaxy spectra are discussed in C03. 

The effect of the age on SSP spectra of the same metal- 
licity (i.e. solar) and IMF (unimodal with slope 1.3) is shown 
in the second row of panels of Fig. 6, where we vary the SSP 
age from 2.0 to 17.8 Gyr. We do not see any significant varia- 
tion either on the Ca n triplet feature or on the overall shape 
of the spectrum. 

4.2 The effects of the IMF 

The third row of panels in Fig. 6 illustrates the effect of the 
IMF on the SSP spectra. In these plots we keep constant the 
age (10 Gyr) and the metallicity (solar) and vary the slope 
of the IMF from [x = 0.3 to 3.3. We see that the third line of 
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Table 3. Predicted CaT* and PaT indices for SSPs of different ages and metallicities (indicated in the first two columns) for different 
IMF shapes and slopes (fi): unimodal, bimodal (according to the definitions given in V96) and the universal and revised IMFs of KOI 
(see § 2.1 and Appendix A). The Salpctcr (1955) solar neighbourhood IMF is given by a unimodal IMF of slope 1.3. The tabulated 
index values were calculated on the basis of the empirical fitting functions given in Paper III (see § 3.3). The spectral resolution is 
1.5A (FWHM). The Paschen contaminated Call triplet index, CaT, can be easily calculated by applying the relation: CaT = CaT* + 
0.93 PaT (see Paper I). 



Unimodal Bimodal K. Universal K. Revised 

[M/H] Age n = 0.3 p = 1.3 \x = 2.3 g = 3.3 ji = 1.3 ji = 2.3 
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the Ca 11 triplet feature is not significantly affected when the 
slope of the IMF varies from p = 0.3 to ~ 2.0, but the line 
rapidly weakens for extremely dwarf-dominated IMFs. This 
result is quantified in Figure 7, where we show the CaT* 
index for different IMF types and slopes. The top panels 
show the values obtained for different slopes and a unimodal 
IMF. The larger the metallicity the larger the weakening of 
the CaT* index as a function of increasing the IMF slope. 



This effect is not remarkable for lower metallicities, although 
we do see an opposite trend for metallicities below [M/H] 
< —1.0. Fig. 4 shows the reason for this behaviour: the larger 
the IMF slope the larger the relative contribution of the MS 
stars (with lower Ca 11 triplet strengths, see Paper III) to the 
total luminosity and to the contribution of the RGB phase 
(see from top to bottom panels) . 

In the two bottom panels of Fig. 7 starting from the 
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Table 3 - continued 



Unimodal Bimodal K. Universal K. Revised 



[M/H] Age /J, = 0.3 n = 1.3 m = 2.3 = 3.3 n = 1.3 ^ = 2.3 
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Figure 6. Model SSP spectra at a = 200 km s -1 and u = 22.2 km s _1 (i.e. the nominal resolution) for SSPs of different metallicitics 
(upper plot), ages (middle) and IMF slopes (bottom). The models were calculated using a unimodal IMF. All the spectra were normalized 
by the flux in the range AA8619-8642A. 



left we show the results for a bimodal IMF for two different 
slopes, i.e. n — 1.3 and 2.3. These plots show that the effect 
of the IMF slope on the CaT* index is much less pronounced 
for this IMF. This result is expected since for the bimodal 
IMF the contribution of the very low MS stars is decreased 
with respect to a unimodal IMF of the same slope. Finally, 



in the last two panels we plot the results for the universal 
and revised IMF of K01. This figure and the index values 
tabulated in Table 3 indicate that the CaT* values obtained 
from the unimodal and bimodal IMFs of slope /i — 1.3 and 
the two IMFs of K01 are very similar. However, the CaT* 
values obtained with the universal IMF match closer the 
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Figure 7. The CaT* index as a function of the IMF. All the models shown here are calculated on the basis of the fitting functions. 
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Figure 8. Broadening correction ACaT*/CaT* and ACaT/CaT 
for the Ca II triplet feature in SSPs. To calculate these corrections 
we have broadened the model spectra by convolving with gaus- 
sians from a = 25 km s _1 up to a = 400 km s _1 in steps of 
25 km s _1 . AI/I is zero for a = 22.2 km s 1 (the spectral res- 
olution of the stellar library). We plot with different line types 
the broadening corrections for a set of representative models of 
different ages (1, 12.59 Gyr), metallicities (0.0, -0.68) and IMFs 
(unimodal, fi = 1.3, fi = 2.8). In grey we show the region cov- 
ered by the broadening corrections obtained for the whole model 
SSP spectral library. The upper envelope for the two indices cor- 
responds to the model /j. = 3.3, [M/H] = +0.2 and 17.78 Gyr, 
whereas the lower envelope refers to the model fi = 0.3, [M/H] 
= -1.3 and 14.12 Gyr. 

ones of the bimodal IMF, whilst those of the revised IMF 
match those of the unimodal shape. This result is easily 
understood by looking at Fig. 1. 



5 CAII TRIPLET FEATURE DEPENDENCE 
ON SPECTRAL RESOLUTION 

In order to study the sensitivity of the Ca II triplet indices to 
the spectral resolution, or galaxy velocity dispersion broad- 
ening (a) , we have broadened the whole SSP model spectral 
library by convolving with gaussians from a — 25 km s _1 up 
to a = 400 km s _1 (in steps of 25 km s _1 ). The Call triplet 
indices were measured for the full set of broadened spectra 
and we fit, for each model, a third-order polynomial to the 
relative changes of the index values as a function of velocity 
dispersion 

1(a) - IN =a + ba + c(J 2 + da s ^ p(g)) (g) 
1(a) 

where ao = 22.2 km s _1 is the nominal resolution of the 
models (FWHM = 1.50 A). Finally, making use of this for- 
mula it is straightforward to relate the indices measured at 
two different velocity dispersions, ai and cri. 

Table 4 tabulates the derived coefficients for the CaT* 
and CaT indices measured over a number of representative 
SSP models. The velocity dispersion corrections for the PaT 
index can be inferred from those of the CaT* and CaT in- 
dices. We plot in Figure 8 the obtained ACaT*/CaT* and 
ACaT/CaT values for this set of representative models. The 
grey region represents the locus of broadening corrections for 



the whole SSP model spectral library. Note that the CaT in- 
dex is significantly more sensitive to the velocity dispersion 
than the CaT* index. This result is in agreement with that 
obtained in Paper I for stars. In that paper, it was shown 
that PaT is quite sensitive to velocity dispersion broadening. 
The improvement in the CaT* sensitivity is explained since 
the effects of the broadening on CaT and PaT indices are 
partially compensated when the CaT* index is computed. 
Paper I also shows that the CaT* index is less sensitive to 
resolution than most of the other popular index definitions 
for this feature. 

The fact that the broadening correction depends on the 
model parameters shows us to what extent one can match 
the effects of galaxy velocity broadening by using a number 
of stellar templates. In fact, the stellar spectra show even 
larger variations than those shown by the models (as it can 
be seen by comparing Fig. 8 to Fig. 5a of Paper I). Note 
that this method has been widely used in the literature to 
match as well the resolution requirements of the Lick/IDS 
system of indices in the optical spectral range in order to 
compare galaxy line-strength measurements with the index 
predictions of models based on that system (Worthey 1994; 
V96). The reader is referred to the extensive review of W097 
for a detailed description of the method. It is evident that 
the method that we propose here to smooth the new model 
spectra to galaxy velocity dispersion (and resolution of the 
data) does not suffer from these broadening correction un- 
certainties. The advantages of this approach over the tradi- 
tional method, based on the stellar templates, has also been 
proven by Falcon-Barroso et al. (2003), who have used the 
new SSP models presented here as templates for deriving 
accurate stellar velocities, velocity dispersions and higher 
order Gauss-Hermite moments profiles. 

It is interesting to see whether the information provided 
by the CaT* index, in terms of the relevant parameters of 
the stellar populations, is varying as a function of galaxy 
velocity dispersion. For this purpose we plot in Figure 9 the 
CaT* index measured on the model SSP spectral library 
smoothed from a = 50 to 300 km s _1 . We see some relative 
variations between the model lines of different metallicities, 
e.g. the separation between the lines corresponding to [M/H] 
= -0.68 and [M/H] = +0.2 is larger at a = 300 km s" 1 
than at a = 50 km s _1 . However, these differences are small 
in comparison to typical observational errors (see Paper I 
and C03) and, therefore, the stellar population parameter 
estimates obtained on the basis of this index do not depend 
significantly on the adopted resolution. 



6 SOLAR METALLICITY MODELS IN THE 
AGE INTERVAL 0.1-1 GYR 

Although the predictions for young stellar populations will 
be described in detail in a forthcoming paper, in this section 
we extend the age of our SSP models down to 0.1 Gyr for 
solar metallicity. When compared to old stellar populations 
here are several major changes in the relative contributions 
to the total luminosity of the different evolutionary stages in 
this age range. In fact, the MS, HB and AGB phases increase 
their contributions, whilst the RGB mostly disappears. In 
Figure 10 we plot a set of representative SSP spectra corre- 
sponding to 0.2 and 2 Gyr for solar metallicity and unimodal 
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Table 4. Coefficients of the broadening correction polynomials (AX/ X = a + bcr + ccr 2 + dcr' 3 ) for the CaT* and CaT indices for a set of 
representative SSP models (//,[M/H],t). 
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Model 


a(xl0~ 3 ) 
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CaT 
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-7.750521 


-4.221853 




1.3, 0.0, 12.6 


-6.121178 
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-3.505068 
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-6.122020 
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Figure 9. The CaT* index measured on the model SSP spectral 
library smoothed to a = 50, 100, 200 and 300 km s _1 . We use a 
unimodal IMF of slope 1.3. 
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IMF of slope 1.3 (for two different spectral resolutions). As 
it was expected, there is a significant strengthening of the 
Paschen series towards younger ages (see for example the 
change in the strength of the lines at either side of the third 
Can triplet line). This result is a consequence of the in- 
creased contribution of hotter MS stars to the total luminos- 
ity in this spectral range. We also see that the continuum 
around the Call triplet feature is significantly affected by 
the molecular band absorptions, which constitute a typical 
feature of the M-type stars. This continuum shape with its 
characteristic slope is attributed to the contribution of the 
AGB phase (see Appendix C) , whilst similar shapes are ob- 
tained for older stellar populations due to the contribution 
of the RGB. 



Figure 10. From top to bottom: SSP model spectra of 0.2 and 
2.0 Gyr for solar metallicity and unimodal IMF with slope 1.3. 
The left panel spectra were smoothed to a = 200 km s _1 , whilst 
the right panel spectra are kept at the nominal resolution of the 
models, i.e. a = 22.2 km s~ 1 . 



In Figure 11 we show the time evolution of the PaT, 
CaT* and CaT indices. The plot shows the strengthening of 
the PaT index and the weakening of the CaT* index towards 
younger SSPs. Interestingly, these changes are compensated 
in a such a way that the CaT index is virtually constant in 
this age range, and have similar values to the ones obtained 
for older stellar populations. Finally, it is worth noting that 
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Table 5. Predicted CaT* and PaT indices for intermediate-age 
SSPs for solar metallicity and different IMF shapes (for the uni- 
modal and bimodal IMFs we adopt //=1.3). The tabulated index 
values are calculated on the basis of the empirical fitting functions 
given in Paper III. The spectral resolution is 1.5A (FWHM). The 
CaT index can be obtained from the relation CaT = CaT* + 
0.93 PaT 



Age 
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CaT* PaT CaT* PaT 



K. Universal 
CaT* PaT 



K. Revised 
CaT* PaT 
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5.056 
4.899 
4.720 
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6.031 



3.952 
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3.796 
3.657 
3.447 
3.202 
3.002 
2.774 
2.531 
2.305 
2.078 



5.755 
5.560 
5.056 
4.898 
4.720 
4.596 
4.502 
4.471 
4.482 
4.473 
4.497 
4.578 
4.656 
4.791 
4.993 
5.198 
5.350 
5.520 
5.703 
5.865 
6.031 



3.953 
4.010 
3.473 
3.586 
3.742 
3.848 
3.963 
4.043 
4.011 
3.996 
3.970 
3.883 
3.797 
3.658 
3.448 
3.203 
3.003 
2.775 
2.531 
2.306 
2.078 



5.755 3.942 
5.561 3.999 
5.057 3.466 
4.901 3.578 
4.723 3.732 
4.600 3.836 
4.507 3.950 
4.477 4.028 
4.488 3.996 
4.479 3.979 
4.504 3.952 
4.585 3.865 
4.662 3.778 
4.797 3.638 
4.998 3.429 
5.202 3.184 
5.353 2.985 
5.522 2.757 
5.704 2.515 
5.865 2.290 
6.029 2.064 



the Paschen series cannot be detected for SSPs of ages larger 
than ~ 1.5 Gyr. In Table 5 we list the predicted CaT* and 
PaT index values corresponding to this age interval for sev- 
eral IMF types. 



7 COMPARISON WITH PREVIOUS 
PREDICTIONS 

In Figure 12 we compare our CaT index predictions (which 
superseed those of V96) with those of MGV (which super- 
seed those of Garcia- Vargas, Molla & Bressan 1998) and 
ITD. These authors employ different Can triplet index def- 
initions (see Paper I for an extensive comparison between 
different systems). However, we can transform their predic- 
tions to our system on the basis of the relations given in 
Table 8 of Paper I. In that Paper we presented two sets 
of calibrations to convert between different systems. The 
first set accounted for the effect of different Can triplet in- 
dex definitions, while the second set can be used to correct 
from flux calibration effects and differences in spectral reso- 
lution. The composite calibrations to convert the predictions 
of MGV (CaT(MGV), which are on the system by Diaz et 
al. (1989), and ITD (CaT(ITD), which employ the system 
by AZ88, to our CaT index are the following: 



oL 
0.1 



1 .0 
Age (Gyr) 



Figure 11. Time evolution of the CaT (dashed line), CaT* 
(solid line) and PaT (dotted line) indices for intermediate aged 
SSPs. These predictions are calculated on the basis of a unimodal 
IMF of slope 1.3 and solar metallicity. The spectral resolution is 
1.5A (FWHM). 




CaT = 0.533 + 1.103 CaT(MGV) 
CaT = -0.155 + 1.183 CaT(ITD). 



(9) 



10 

Age (Gyr) 



Figure 12. Comparison of the CaT index predictions by different 
authors for SSPs of different metallicities and ages. All the models 
are calculated with a Salpeter IMF. The dotted lines refer to the 
predictions of MGV, the dashed lines refer to ITD and the solid 
line represents our predictions. To perform this comparison we 
have transformed their predictions to our system on the basis of 
the formulae given in Paper I (see the text). 
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In Fig. 12 we compare the tabulated predictions of these 
authors (converted to our system) with ours. Among the 
main differences, we must remark that their CaT index pre- 
dictions increase as a function of metallicity for all metallic- 
ity regimes, whilst we obtain a saturation of this index for 
metallicities above ~ —0.5. This saturation can be explained 
in part by the fact that the larger the metallicity the lower 
is the difference between the contributions of the Main Se- 
quence (lower CaT values) and the Red Giant Branch (larger 
CaT values) to the total luminosity of an SSP, as it can be 
seen in Fig. 4 (second row of panels, i.e. fi = 1.3). Moreover, 
in metal-rich SSPs the RGB is populated with an increas- 
ing fraction of cold stars, and, for temperatures cooler than 
~ 3750 K, the CaT index decreases dramatically with de- 
creasing temperature (see Fig. 7 of Paper III). The main 
point is that this behaviour for the M giants is not pre- 
dicted by the fitting functions used by these authors since 
their stellar libraries do not include these stars. 

In the case of ITD, the coolest stars included in their 
library are of K3 spectral type. MGV extrapolate to lower 
temperatures the theoretical fitting functions of j0rgensen et 
al. (1992), valid only for T off > 4000 K. Note that although 
MGV try to justify this extrapolation in their Figure 1, their 
argument is not completely correct. For instance, for a giant 
star of T c ff = 3000 K (which has a gravity of logg ~ 0), 
the J0rgensen et al. extrapolated fitting functions predict 
CaT = 9.84 A (11.39 A when corrected to our system), 
whereas the actual CaT measured values in our library stars 
of similar parameters range from 2.3 to 5.0 A. 

These effects may also explain the fact that these au- 
thors obtain slightly larger CaT values for increasing ages, 
whilst our predictions for the time evolution of the CaT in- 
dex is towards slightly lower values for SSPs of large metal- 
licities. 

For lower metallicities, our CaT predictions are in bet- 
ter agreement with the predictions of ITD, whilst those of 
MGV show much lower values. This can be explained on 
the basis of Figure 12 of Paper III, where it can be seen 
that the fitting functions of J0rgensen et al. predict a larger 
dependence on metallicity than what is actually observed. 

The first paper in synthesizing SSP spectra at high res- 
olution in this spectral range is that of SBB. These authors 
follow a completely different approach than the one used 
in this paper. They employ a fully synthetic stellar spec- 
tral library calculated by Schiavon & Barbuy (1999), which 
is based on model photospheres and molecular and atomic 
line lists. Another important difference is that they adopt 
the old Padova set of isochrones (B94) , whilst we make use 
of the new Padova set. Figure 13 shows a comparison of a 
representative SSP model spectrum, kindly provided to us 
by these authors, to an equivalent spectrum of our model 
library. These spectra correspond to a Salpeter IMF, solar 
metallicity and 13 Gyr and have a resolution 2 A (FWHM). 
They were normalized according to the flux in the spectral 
range covered by our models. In the lower panel of Fig. 13 we 
plot the difference between these two models. Although we 
see some differences in the pseudocontinua, particularly in 
the bluest one, the main difference is in the depth of the Ca II 
triplet feature. We have measured the CaT index of the SBB 
spectra and find that their Ca II triplet strengths are 2.5-3 A 
lower than the values we obtain for the SSP spectra of our 
model library. In particular, for the spectra plotted in Fig. 13 




8500 8600 8700 

Wavelength (A) 



Schiovon et ol (2000) 
This paper 



8800 8900 




8500 8600 8700 



8800 8900 



Figure 13. Normalized representative SSP model spectra of SBB 
(thin line) and this work (thick line). The two models have a 
Salpeter IMF, solar metallicity and 13 Gyr. The spectra were 
smoothed to FWHM = 2A. The difference between these two 
spectra is shown in the lower panel. 



we obtain 4.25 A for the SBB spectrum and 7.33 A for our 
spectrum. Schiavon (private communication) suggests that 
the difference in the Ca II triplet strength can be attributed 
in part to the fact that their synthetic stellar spectral library 
was calculated without adopting a NLTE approach, which it 
is needed to properly model the depth of the lines. In fact we 
have compared the empirical stellar spectra of our library to 
synthetic stellar spectra of similar atmospheric parameters 
of their model grid and found that their Ca II triplet values 
are systematically lower for the giant stars. These differences 
may be as high as 5 A. On the other hand, we find a reason- 
ably good agreement for dwarf stars. This comparison shows 
that further work on the theoretical libraries, including ap- 
propriate input physics and updated opacities, is required 
for this spectral range in order to be able to make use of 
these libraries for stellar population synthesis modeling. 



8 COMPARISON TO GLOBULAR CLUSTERS 
AND EARLY- TYPE GALAXIES 

8.1 Globular clusters 

Globular clusters are the ideal candidates to check our model 
predictions since they can be treated as single-age, single- 
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Figure 14. Galactic globular cluster Can triplet index measure- 
ments given by AZ88 (i.e. ^ W) versus the metallicities predicted 
by our models, when transforming from their system to ours, by 
means of the CaT index. We assumed a Salpeter IMF and that 
the clusters are 14 Gyr old (see the text). For clusters with [Fe/H] 
< —1.68 we extrapolated our model predictions. The solid line 
represents the AZ88 calibration. 



metallicity, stellar populations. We follow here the test sug- 
gested by ITD to calibrate our Ca II triplet predictions ver- 
sus the semi-empirical metallicity calibration of AZ88. For 
this purpose, we need to transform their galactic globular 
cluster Call triplet index measurements (i.e. ^ W , follow- 
ing the notation of these authors) to our system by means 
of the CaT index. This conversion is required since these au- 
thors used different spectral resolution and dispersion and 
did not correct their data to a relative flux scale. Therefore 
we first applied the broadening correction given in Table 7 
of Paper I to translate their measurements, i.e. at resolution 
4.8 A, to ours (1.5 A). The second step is to transform their 
W index (on their instrumental response curve) to our 
CaT index (on a relative flux scale) by making use of Eq. 9. 

The fact that our models show that the Call triplet 
feature exhibits a negligible dependence on the age of the 
stellar population, for the age range that is usually assumed 
for the globular clusters, makes it possible to neglect this pa- 
rameter. Figure 14 shows our metallicity estimates obtained 
from a comparison of the corrected CaT measurements to 
our predictions for models of 14 Gyr and unimodal IMF 
of slope 1.3. We note that we extrapolate the models for 
[Fe/H] < -1.68. The solid line represents the AZ88 cali- 
bration ([Fe/H] CaII = -4.146 + 0.561 X)W). Fig. 14 shows 
the good agreement achieved for metallicities smaller than 
~ —0.5. This result probes the validity of using the Call 
triplet of integrated spectra as an alternative method to de- 
rive metallicities in globular clusters of [M/H]< —0.5. 

Fig. 14 also shows a small shift between the obtained 



metallicities A [Fe/H] < 0.15 dex. In part, this shift may 
be explained by the large sensitivity of the ^ W index of 
AZ88 to the spectral resolution, as it was shown in Paper I. 
In fact, we have tested that for a resolution mismatch of 
~ 20 km s _1 the obtained offset vanishes. This mismatch 
could either be originated if the AZ88 resolution differs from 
4.8 A (th ese authors obtained their observations in differ- 
ent telescopes and used different instrumental setups) or by 
our assumption of taking the broadening correction repre- 
senting a M0 star, according to Paper I. An important fac- 
tor to be taken into account for explaining the obtained 
metallicity shift is the fact that we adopted the Carretta 
and Gratton (1997) scale for the cluster stars of our stellar 
spectral library (see Paper II). Moreover, for the stars of 
M 71 we adopted here an even lower metallicity value (i.e. 
[Fe/H] = —0.84) than that predicted by Carretta & Gratton 
(1997) (i.e. [Fe/H] = -0.70), as discussed in § 3.1.1. It has 
been largely discussed in the literature that this scale, based 
on high-dispersion spectra, shows significant deviations from 
the Zinn & West (1984) scale, which was adopted by AZ88. 
This effect is more pronounced for the largest metallicities. 
We refer the reader to the discussion presented in Rutledge 
et al. (1997). 

Finally, we note that for metallicities larger than ~ —0.5 
the measured CaT index does not increase as a function of 
metallicity, departing from the metallicity scale relation of 
AZ88. Therefore, these measurements are in agreement with 
our model predictions. However, as a result of this satura- 
tion, we are not able to provide accurate metallicity esti- 
mates for the most metal-rich globular clusters on the basis 
of the Ca II triplet . 

8.2 Early- type galaxies 

In C03 we have applied these models to a large sample of 
early- type galaxies of different luminosities. In that paper 
we analyse the Can triplet as well as newly defined fea- 
tures, such as the slope of the continuum around the Can 
triplet lines and the Mgi feature at 8807A. We therefore re- 
fer the reader to that paper for a complete description of the 
galaxy sample and for an extensive discussion of their stellar 
populations. In this section, we aim at showing the poten- 
tial use of the new models by selecting two representative 
galaxies from C03 sample that have already been analysed 
in the blue spectral range (V01A) on the basis of very high 
quality spectra that allowed us to measure the new age in- 
dicator H7 CT of VA99, which does not depend on metallicity. 
These galaxies are NGC 4478 and NGC 4365. The galaxy 
spectra that we use here correspond to the central 4 arcsec 
aperture. 

The synthesized model spectra both in the blue (V99) 
and in the spectral region of the Call triplet can be used 
to analyse the observed galaxy spectra in a very easy and 
flexible way, allowing us to adapt the theoretical predic- 
tions to the characteristics of the data, instead of proceed- 
ing in the opposite direction as, for example, we must do 
when we have to compare our data with models based on 
the widely used Lick/IDS system (see V99 for an extensive 
discussion of the advantages of the new approach). After 
smoothing the synthetic SSP spectra, with flux-calibrated 
spectral response curves, to the measured resolution (i.e. 

^total ^galaxy ~t~ ^"instr 

), we can analyse the galaxy spectrum 
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in its own system. The spectral resolution of the near-IR 
spectra that we use here is <7i nB tr = 42 km s _1 . We then mea- 
sure our favorite spectral features, such as the Call triplet, 
in both the models and the data in order to build up index- 
index diagrams that can help us to obtain the most relevant 
stellar population parameters. 

In Figure 15 we combine our new Call triplet index 
measurements (C03) with the blue index measurements of 
V01A. The CaT* and CaT indices are plotted versus the 
H7ioo«r<i75 and H7 2 25«r<300 indices of VA99. Each R^ a 
index definition provides stable and sensitive age predic- 
tions within the a ranges quoted in the subindices. There- 
fore, whereas the first index is the one appropriate to study 
NGC 4478 (i.e. <7 t otai = 140 km s _1 ), the second is adequate 
to study NGC 4365 (i.e. cr tota i = 260 km s _1 ). The metallic- 
ity of the models increases from the left to the right panels 
as quoted above the top panels. We use models of unimodal 
IMF characterized by its slope /i, which is constant along 
the vertical solid lines. We vary this slope, indicated by the 
strength of the line, from /i = 0.3 (thinnest) to /i = 3.3 
(thickest) in steps of A/i = 0.5. The third vertical line start- 
ing from the thinnest (i.e. right) corresponds to the Salpeter 
(1955) (i.e. /i = 1.3). Models of constant age are shown by 
thin dotted horizontal lines. 

All the model grids of Fig. 15 look approximately or- 
thogonal, which indicates that for a given metallicity we 
are able to separate the effects of the age from the ones of 
the IMF slope. According to these plots, smaller U'ja values 
mean larger ages and smaller CaT* or CaT indices mean 
steeper IMF slopes (i.e. dwarf-dominated). We see that the 
inferred ages do not depend on the metallicity of the model 
grid in use due to the insensitivity of H7 CT to this parameter. 
Moreover, the fact that the Can triplet indices are insensi- 
tive to the metallicity for the range of metallicities covered 
in these plots (-0.38 < [M/H] < +0.2) makes it possible 
to provide almost unique IMF slope solutions. All the dia- 
grams indicate a mean luminosity weighted age of ~9 Gyr 
for NGC 4478 and ~ 15 Gyr for the giant elliptical galaxy 
NGC 4365. It is worth noting that the signal-to-noise of the 
NGC 4365 spectrum was not high enough for an accurate 
measurement of the H7 CT index. However, the obtained age 
is not different from the one derived on the basis of the H/3 
index (see V01A), and the age estimate given in Davies et 
al. (2001) from integral field spectroscopy. From the H7 CT 
versus CaT* diagrams we find n ~ 2.3 for NGC 4478 and 
/i ~ 2.8 for NGC 4365, which indicate rather steep IMFs in 
comparison to the Salpeter value. However, it is worth recall- 
ing that variations in the theoretical prescriptions adopted 
for building-up the stellar tracks and isochrones might drive 
to important changes in the predicted IMF slope. In fact, 
if we replace the model grids of Fig. 15 with another set 
synthesized on the basis of the old Padova isochrones (i.e. 
B94), which yield smaller CaT* strengths (i.e. ~ —0.5 A) as 
a result of a slightly cooler RGB phase (see Appendix C2 
for details on the differences in the adopted theoretical pre- 
scriptions), we would have obtained /i ~ 1.3 for NGC 4478 
and n ~ 1.8 for NGC 4365 (the H7 CT age indicator does not 
vary significantly as a result of this replacement, see V01B). 
Therefore these results must be taken into account on a rel- 
ative basis, whereas the discussion about galaxy trends is 
more secure (C03). 

In Figure 16 we plot the CaT* index versus different 



metallicity indicators in the blue spectral range, i.e. [MgFe] 
(Gonzalez 1993), Mgb (W94), Fe3 (K00) and Ca4227 (W94). 
In order to be able to show on a single plot the two galaxies 
we have smoothed the spectrum of NGC 4478 to match the 
resolution of NGC 4365. Each panel shows an approximately 
orthogonal model grid, with the position of the grid varying 
from the left to the right panels as a function of increas- 
ing metallicity. We, therefore, are unable to obtain a robust 
age determination on the basis of this figure. However, if 
we adopt for these galaxies the ages inferred from Fig. 15 
we are in position to obtain the metallicity. The [MgFe] - 
CaT* diagrams suggest [M/H] [MgFe] ~ +0.1 for NGC 4478 
and solar metallicity for NGC 4365. From the Mg& - CaT* 
diagrams we obtain [M/H]M g b ~ +0.2 for the two galaxies. 
The Fe3 - CaT* plots suggest [M/H] Fo3 ~ 0.0 for NGC 4478 

and [M/H] F c3 0.1 for NGC 4365. Finally, from the 

Ca 4227 - CaT* diagrams we obtain [M/H] Ca 4227 ~ -0.2 for 
NGC 4478 and [M/H] C a4227 ~ -0.4 for NGC 4365. Interest- 
ingly, despite the fact that NGC 4365 is a larger galaxy, its 
metallicity seems to be slightly lower than the one obtained 
for NGC 4478. 

Fig. 16 shows us that the obtained metallicities are 
strongly influenced by non-solar abundance ratios in agree- 
ment with previous determinations (e.g. see, for recent refer- 
ences, K00; TOO; V01A, and PS02). Particularly interesting 
is the fact that the lowest metallicities are inferred when us- 
ing the Ca 4227 - CaT* diagrams. Despite the fact that Ca, 
like Mg, is an a-element, Ca 4227 does not track Mg6 (see 
V97; Peletier et al. 1999; TOO; V01A; PS02). Although this 
result is in disagreement with our present knowledge of the 
nucleosynthesis theory (Woosley & Weaver 1995), if these 
galaxies were deficient in Ca, causing the Ca 4227 line to 
be smaller than predicted by scaled-solar models, we should 
consider the possibility that the Can triplet feature might 
be reflecting this deficiency as well. This suggests, as an al- 
ternative scenario, that the IMF is Salpeter-like and the low 
Can triplet values are due to Ca deficiencies rather than 
from the steepening of the IMF. 

A possible approach that can be followed for interpret- 
ing these results is to build-up SSP predictions based on 
non solar abundance ratios. Such isochrones were predicted 
by, e.g. SW98, S00, VandenBerg et al. (2000), Kim et al. 
(2002), for several a-enhancement ratios. These calculations 
were motivated by the finding that elliptical galaxies show 
an enhancement of Mg over Fe when compared to scaled- 
solar stellar population model predictions (e.g. Peletier 1989; 
Worthey et al. 1992; K00; V01A). However to predict stellar 
population models of different a-enhancement ratios we re- 
quire stellar tracks built-up on the basis of appropriate opac- 
ity tables and energy generations, and the corresponding a- 
enhanced stellar spectral libraries. It is not yet clear what 
ratios should be adopted for the different a-elements. In fact 
such a-enhancement stellar models have been calculated by 
various authors assuming a constant enhancement for each 
a-element (e.g. VandenBerg et al. 2000; Kim et al. 2002), 
or somewhat more empirically-based element mixture (e.g. 
SW98; S00). Furthermore Ca 4227 does not track Mg (e.g. 
V97). However, such theoretical stellar spectral libraries are 
not yet available and all the empirical libraries such as the 
one used here mostly follow the Galactic disk element ra- 
tios, particularly for the metallicity regime characteristic of 
the galaxies shown in Fig. 16 (i.e. around solar). Amongst 
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Figure 15. Plots of CaT* and CaT indices versus the H-y a age indicators of VA99. Filled symbols represent the index measurements 
for the central apertures (4 arcsec) of two galaxies in common between the samples of C03 (for the Can triplet indices) and V01A (for 
the H7 indices). The triangle represents NGC 4478, whilst the square represents NGC 4365. Overplotted are the models by V99 for the 
H7 indices and those from this paper for the Call triplet feature. All index measurements were performed after smoothing the model 
SEDs to the appropriate galaxy velocity dispersion and instrumental resolution. We increase the metallicity of the model grids from left 
to right panels as quoted in the upper plots. We use a unimodal IMF, varying in each panel its slope from fi = 0.3 (thinnest vertical 
line) to fi = 3.3 (thickest) by steps of Afj, = 0.5. The third vertical line starting from the thinnest (i.e. right) corresponds to the Salpeter 
(1955) value (i.e. fj, = 1.3). Finally, thin dotted horizontal lines represent models of constant ages, which are quoted in gigayears. 
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Figure 16. Plots of CaT* versus different metallicity indicators in the blue spectral region. Lines and symbols have the same meaning 
as in Fig. 15. In these plots all index measurements were obtained at resolution <r tota i = 260 km s _1 . The blue indices are taken from 
V01A. 
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the most important effects of adopting such isochrones is 
the fact that, for a given total metallicity, the a-enhanced 
mixtures lead to lower opacities, which translates into an 
increase of the temperature of the stars in both the Main 
Sequence (MS) and the RGB phases. However, for a given 
[Fe/H] metallicity, we obtain the opposite trend since the to- 
tal metallicity of the a-enhanced models is larger. Another 
important result from these studies is that, for metal-rich 
stellar populations, a-enhanced isochrones cannot be mim- 
icked with a scaled-solar isochrone of different metallicity. 

Although it would be interesting to have models with 
varying [Ca/Fe] ratios as well, stellar libraries and isochrones 
calculated with such abundance ratios are not yet available. 
Moreover, it is not clear how we expect [Ca/Fe] to behave. It 
is worth noting that we did not obtained any significant cor- 
relation of the [Ca/Fe] abundance ratios with the residuals 
of the predictions of the fitting functions and the observed 
values for the Call triplet in our stellar library (Paper III). 
Despite the fact that we are not in position to build-up such 
fully self-consistent non-solar element ratios SSP models, we 
have tested in Appendix C2 the effects of adopting the a- 
enhanced isochrones of S00. These isochrones are calculated 
on the basis of the same input physics as that used in GOO for 
old stellar populations. For total metallicity [M/H] > we 
obtain larger CaT* strengths (ACaT* ~ 0.5 A). If we plot 
these model predictions on Fig. 15 we obtain even steeper 
IMF slopes. Alternatively, if a Salpeter-like IMF is used, the 
a-enhanced isochrones lead to a larger calcium underabun- 
dance problem. 

It is however interesting to see how the model spectra 
of the adopted IMF slopes, ages and metallicities match the 
full spectral region around the Ca II triplet feature for these 
galaxies. In Figure 17 we show the spectra of NGC 4478 
(top panel) and NGC 4365 (bottom panel). For NGC 4478 
we overplot several models of age 9 Gyr as suggested by 
Fig. 15. In order to illustrate the effect of the IMF slope we 
have chosen two model spectra corresponding to /u = 1.3 (as 
a reference) and fx — 2.3 (the most suitable fit, according 
to Fig. 15 and Fig. 16). The selected metallicity is around 
solar according to the [MgFe] and Fe3 indices. All the spec- 
tra were normalized according to the flux measured in the 
spectral region AA8619-8642 A. It is apparent that the ref- 
erence model with (j, = 1.3 provides slightly deeper Can 
triplet lines than observed in the galaxy spectrum. We also 
see that the model with fi — 2.3 provides a better fit to the 
whole spectral region, including the slope of the continuum 
around the Call triplet feature. Finally we overplot another 
model spectrum of n = 2.3 and [M/H] = —0.1, which seems 
to provide a somewhat better fit to the overall galaxy spec- 
trum. 

In the lower panel of Fig. 17 we plot the spectrum of 
NGC 4365 and overplot several model spectra of similar 
age (15.85 Gyr) and IMF slope (p = 2.8). The four over- 
plotted models were selected to have the metallicities in- 
ferred on the basis of the four metallicity indicators studied 
in Fig. 16. Overall, we see that the best fits are achieved 
for the metallicities in the range —0.2 < [M/H] < +0.1, in 
good agreement with the predictions obtained on the basis 
of the [MgFe] (i.e. ~ 0.0) and Fe3 (i.e. ~ -0.1) indices. On 
the other hand, the models of metallicities obtained from the 
Mgb (i.e. ~ +0.2) and Ca4227 (i.e. ~ -0.4) provides slightly 
steep and rather flat continua around the Call triplet fear 
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Figure 17. Spectra of NGC 4478 (upper panel) and NGC 4365 
(lower panel). Overplotted are several model spectra selected on 
the basis of Fig. 15 and Fig. 16 (see text for details). All the 
spectra were normalized according to the flux measured in the 
spectral region AA8619-8642A 



ture respectively. This result suggests that the metallicity 
inferred on the basis of the Fe3 and [MgFe] indices provide 
reasonably good fits to the overall spectrum in this spectral 
range. 



9 CONCLUSIONS 

We have presented a new evolutionary stellar popula- 
tion synthesis model in the near-IR spectral region cover- 
ing the range AA8350-8950 A at a spectral resolution of 
1.5 A (FWHM). The ultimate aim of the models is to derive 
reliable predictions of the Call triplet strength for stellar 
populations in a wide range of ages, metallicities and IMF 
shapes. Apart from spectral energy distributions, the model 
predicts measurements in a set of new line indices, namely 
CaT*, CaT and PaT, which we have defined in Paper I to 
overcome some of the limitations of previous Call triplet 
index definitions in this spectral range. In particular, the 
CaT* index is specially suited to remove the contamination 
from Paschen lines in the integrated spectra of galaxies. 
The stellar population model presented here is a revised 
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version of the model of V96 and V99. Several aspects of the 
model have been updated in order to produce a state-of-the- 
art output. The isochrones of B94 have been replaced by 
the newest Padova isochrones (GOO), based on solar abun- 
dance ratios. The theoretical parameters of the predicted 
stars have been transformed to fluxes and colours on the 
basis of almost fully empirical relations based on extensive 
photometric stellar libraries. Apart of the IMF shapes used 
in V96 we use the two multi-part power-law IMFs recently 
proposed by KOI. The main ingredient for the predictions 
presented here, is the new extensive empirical stellar spectral 
library presented in Paper I and II, from which a subsample 
composed of more than 600 stars has been carefully selected. 
The sample shows an unprecedented coverage of the stellar 
atmospheric parameters (T e tt, log g and [M/H]) for stellar 
population synthesis modeling. 

Two main products of interest for stellar population 
analysis are presented: i) a spectral library for SSPs with 
metallicities —1.7 < [Fe/H] < +0.2, a large range of ages 
(0.1-18 Gyr) and IMF types, and it) line-strengths calculated 
on the basis of the empirical fitting functions presented in 
Paper III for the CaT* , CaT and PaT indices, and for which 
conversion formulae to previous definitions can be found in 
Paper I. Tables and model spectra are available electroni- 
cally at the web pages given in § 3, together with the Can 
triplet stellar library. 

The newly synthesized model spectra can be used to an- 
alyze observed galaxy spectrum in a safe and flexible way, 
allowing us to adapt the theoretical predictions to the na- 
ture of the data instead of proceeding in the opposite direc- 
tion. The synthetic SSP library, with flux-calibrated spec- 
tral response, can be smoothed to the same resolution of 
the observations or to the measured internal galaxy velocity 
dispersion, making it possible to fully use all the informa- 
tion in the data. This opens the way (as do the models of 
V99 for the blue spectral region) for new applications. For 
example, if part of the spectrum is corrupted or affected by 
strong sky-lines, one can easily define new absorption lines 
on both models and data. Also, these models can be used 
as templates for determining stellar kinematics in galaxies 
(see Falcon-Barroso et al. 2003). Moreover, the models offer 
us a great opportunity to accurately study galaxies at larger 
redshifts. 

We have analyzed in detail the behaviour of the Can 
triplet feature in old-aged SSPs, finding the following re- 
sults: i) the strength of the CaT* index does not vary much 
for ages larger than ~ 3 Gyr, for all metallicities, ii) this 
index shows a strong dependence as function of metallicity 
for values below [M/H] ~ —0.5 and Hi) for higher metallic- 
ities, this index does not show any significant dependence 
either on age or on the metallicity, being more sensitive to 
changes in the slope of power-like IMF shapes. These models 
open up the analysis of a suitable dwarf-giant discriminator 
in the analysis of stellar populations in galaxies. It is worth 
noting that the saturation of the Call triplet feature for 
metallicities above —0.5 has not been predicted by previous 
models. This prediction might be supported by the most 
metal-rich Galactic globular clusters and galaxy data ana- 
lyzed here. Furthermore, it might also be supported by the 
measurements of the CaT* and CaT indices in large sam- 
ples of elliptical galaxies (Saglia et al. 2002; C03) and bulges 
(Falcon-Barroso et al. 2002). 



For stellar populations in the age interval 0.1-1 Gyr we 
find that the Paschen lines, as defined by the PaT index, 
become more prominent and the CaT* index decreases sig- 
nificantly. However the two effects are compensated in the 
CaT in such a way that that this index is virtually constant 
along this age range. Another interesting result is that the 
overall shape of the continuum around the Ca II triplet fea- 
ture for large metallicities resembles more that of M-type 
stars rather than that of K-giants (as it is the case for the 
blue spectral range). 

By comparing globular cluster data (§ 8.1) we show that 
the derived metallicities are in excellent agreement with the 
ones from the literature, probing the validity of using the 
integrated Call triplet feature for determining metallicities 
for these stellar systems. It is worth noting that this is possi- 
ble since the models predict that the Ca II triplet indices are 
virtually constant for stellar populations of ages larger than 
~3 Gyr and strongly dependent on metallicity for [M/H] 
< -0.5. 

In § 8.2 we have applied the models to two early-type 
galaxies of different luminosities, NGC 4478 and NGC 4365 
(from C03), for which a detailed stellar population analysis 
based on the optical spectral range was published in V01A. 
We propose several index-index diagrams, i.e. CaT* versus 
several age indicators, and CaT* versus several metallicity 
indicators, which provide virtually orthogonal model grids 
from which to disentangle, unambiguously, the relevant pa- 
rameters of their stellar populations. The Ca-abundance as 
inferred from the Call triplet does not follow Mg, as nucle- 
osynthesis calculations predict, in agreement with measure- 
ments of the Ca 4227 line in the blue (V97). We should, 
however, note that Ca II triplet measurements cannot be fit- 
ted unless a very dwarf-dominated IMF is imposed, or unless 
the Ca-abundance, as measured from the Can triplet, is even 
lower than the Fe-abundance. However, it is worth recalling 
that variations in the theoretical prescriptions adopted for 
building-up the stellar tracks and isochrones might lead to 
significantly lower IMF slope estimates and, therefore, these 
results must be taken into account on a relative basis, be- 
ing the discussion about galaxy trends more secure (see e.g. 
Saglia et al. 2002; C03). Moreover if we adopted a-enhanced 
isochrones we would obtain steeper IMF slopes compared to 
scaled-solar. Alternatively, if a Salpeter-like IMF is used, the 
a-enhanced isochrones lead to a larger calcium underabun- 
dance problem. It is also important to note that the overall 
shape of the spectrum around the Ca II triplet feature and in 
particular the characteristic slope of the continuum is signif- 
icantly better represented by choosing SSPs of metallicities 
as inferred from Fe or [MgFe] indices, rather than the ones 
based on the Mgb or Ca 4227 lines. More details can be found 
in C03, where a large sample of ellipticals are discussed in 
detail revealing very interesting trends. 
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APPENDIX A: THE IMF TYPES 

In this appendix we summarize all the IMFs used in this 
paper, which include the two IMF shapes given by V96, i.e. 
unimodal and bimodal, and the new IMFs proposed by K01, 
i.e. universal and revised. These four IMF shapes are plotted 
in Fig. 1. 

The unimodal IMF has a power-law form characterized 
by its slope /j as a free parameter 

:/3m- ( " +1) . 



<f>(m) 



(Al) 



Therefore, the Salpeter (1955) solar neighbourhood IMF is 
obtained when /i = 1.3. 

The bimodal IMF is similar to the unimodal case for 
stars with masses above O.6M , but decreasing the weight 
of the stars with lower masses by means of a transition to 



a shallower slope, which becomes flat in the log($(log(m))) 
- log(m) diagram for masses lower than O.2M0. The main 
motivation for this IMF is to achieve a reasonable good fit 
to the observational data of Scalo (1986) and Kroupa, Tout 
& Gilmore (1993) on the basis of a single free parameter, 
i.e. fi, rather than using several segments as it is the case for 
the IMFs of K01. The main advantage of the bimodal IMF 
is that it allows us to vary the parameter \i in the same way 
as for the unimodal case. The form of this IMF^ is given by 



$(m) = (3 < p^m) 



rn < mo 

mo < m < m 2 



(A2) 



m (m+1) , m > m 2 , 



where mo, mi and mi are 0.2, 0.4 and 0.6 Mq respectively. 
p(m) is a spline for which we calculate the corresponding 
coefficients solving for the following system obtained by the 
boundary conditions 



p(m ) 
p'(m ) 

p(m 2 ) 
p'(m 2 ) 



(A3) 
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We have introduced in this paper the two IMF shapes 
recently proposed by K01. The universal IMF is a multi-part 
power-law IMF, which has the following form 



$(m) 




, m < mo 

, mo < m < mi 

, m > mi, 



(A4) 



where mo and mi are 0.08 and 0.5 M©, respectively. 

We have also implemented the revised multi-part power- 
law IMF of K01, which tries to correct for the systematic 
effects due to unresolved binaries on the single-star IMF. 
The main effect is that the slope of the IMF is steeper 
than the universal IMF by A/i ~ 0.5 for the mass range 
0.08 < m < 1 M Q . Its shape is summarized in the following 
equation: 



$(m) = P < 



(sj) 18 ,m <m<mi 
(^r^fer 2 ' 7 ,mi<m<m 2 



where mo, mi and m 2 are 0.08, 0.5 and 1 Mq, respectively. 

Finally the constant P is calculated via the normaliza- 
tion 



r 

P / <E>(m)dm = 1, 

J mi 



(A6) 



where mi and m u are the lower and upper mass cutoffs. 
We adopt 0.01 M Q and 120 M Q for these stellar masses, 
respectively. 



t Note the erratum present in Eq. 4 (case rn < O.2M0) of V96, 
which should be divided by m. 
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APPENDIX B: COMPUTATION OF A 
REPRESENTATIVE STELLAR SPECTRUM 

We provide here a detailed description of how we interpolate 
spectra in the database and compute the spectrum corre- 
sponding to a given set of requested parameters (8 , log go 
and [M/H]o). The most straightforward approach to follow is 
to let the code find all the stars enclosed within a given cube 
centered on these atmospheric parameters. However, since in 
many cases there is a lack of symmetry in the distribution 
of stars around a given point in the 3-parameter space, we 
have divided the original box in 8 cubes all with one corner 
at do, log go, [M/H]o. In this way we make sure that there 
are stars in all the directions with respect to the required 
point. The code finds all the stars Nj enclosed within a cube 
j (where j = 1, . . . , 8) and combines their spectra according 
to 



J,([M/H] 



(Bl) 



where W/ is the weight assigned to the star i within the 
cube j: 
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where SN Ui is a measure of the quality of the empirical 
spectrum Sx? estimated by its signal-to-noise normalized in 
the following way 
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where SNgg.7 is the limiting SN value where the cumula- 
tive SN distribution fraction of our stellar sample reaches 
the upper 99.7 percentile (]T SN^SNi < SNg 9 . 7 )/ £ SNi). 
We note, however, that a spectroscopic binary or a star with 
anomalous signature of high variability that is present in the 
stellar library (see § 3.1.1) would be able to contribute sig- 
nificantly within a given cube if its SN Ui is high. To prevent 
this possibility we have chosen to get down the SN ni of these 
stars to the value where the cumulative SN distribution frac- 
tion of our stellar sample reaches the lower 0.3 percentile. 

The term within the square brackets in Eq. B2 repre- 
sents the weight resulting from the position of the star i in 
the parameter space (8, log 3, [M/H]). We have assumed a 
gaussian-like function, which assigns larger weights to the 
stars closer to the requested point (9o, log go, [M/H]o). This 
term, as well as SN ni , varies from to 1. The ag , o"io g90 
and <7[m/h] values (generically u P0 ) are estimated on the 
basis of the density of stars at the requested point in the 
parameter space po- We assumed an inversely proportional 
gaussian-like function of the form: 



and 



1 I PQ-PM V 
, 2 V CT P0 / 



(B4) 



where pm is the maximum density. Assuming that the small- 
est a Pln is obtained at the point of maximum density (pm) 
and that the largest a PM is reached when p — > we derive 
a PQ and then Eq. B4 can be written as: 



Figure Bl. Grey levels representing the values of ctq, <Xi gg 
CT [M/H] (from top to bottom plots), calculated for [M/H] = —0.7 
(left) and [M/H] = 0.0 (right) in the 9 — \ogg plane. Black means 



<r e = 0.009, <r lo 



0.18 and <T[m/h] = 0.09, whereas white 



means erg = 0.17, o"i og9 = 0.51 and crj M/ / H ] 
sponding plots. 
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We have assumed <7e m , f"i og9m and cr[M/H] m to be the 
minimum uncertainty in the determination of 8, logg and 
[M/H] respectively. We adopted a 0ln = 0.009, o"i og g m = 0.18 
and cr[M/H] m = 0.09 according to the values given in Pa- 
per II. It is worth noting that the points where the uncer- 
tainties are the smallest coincide with the most populated 
regions of the parameter space, e.g. dwarfs of ~ 5800 K 
and giants of ~ 4800 K of solar metallicity. On the other 
hand, for a$ M , ciogg M and i7[m/h] m we assumed a range of 
1/10 of the total covered by these parameters in our stel- 
lar library. We adopted ctom = 0.17, o"io g9M = 0.51 and 
C[m/h] m = 0.41. For og we also assumed the condition that 
60 < T of f < 3350 K (corresponding to the adopted a e limit- 
ing values when transforming to the T e ft scale). 
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At a given point po is calculated by counting the number 
of stars present in a box of size 3xu Pm . If no stars are found, 
the box is simultaneously enlarged in its three dimensions, in 
steps of lxcTp m , until at least one single star is reached. For 
stars with temperatures larger than 9000 K and lower than 
4000 K we do not take into account the metallicity parame- 
ter, since its uncertainty is very large and since we adopted 
[Fe/H] = for stars with unknown metallicity in these tem- 
perature ranges. In practise, we adopt pM as the density 
value where the cumulative density fraction reaches the up- 
per percentile 99.7. For larger densities we adopt this value 
as well. An example of the resulting ag, a\ ogg and a [m/h] is 
shown in Figure Bl for two values of the metallicity. Black 
means ag — 0.009, o"i og9 = 0.18 and <7[m/h] — 0.09, whereas 
white means ag — 0.17, a"i og9 = 0.51 and <7[m/h] = 0.41 
in the corresponding plots. Overall, the grey levels for the 
plots corresponding to [M/H] = 0.0 are significantly darker 
than the ones for [M/H] = —0.7 due to the fact that our 
library is more populated with solar metallicity stars. Note 
that the MS and RGB phases are clearly visible. An inter- 
esting feature is that for the ag plots the grey levels become 
gradually darker towards hotter stars as a result of the fact 
that we assumed that the largest ag value cannot be larger 
than 1/10 of total range in the T e s scale. Note the jump in 
the grey levels when 6 reaches values around 0.55 and 1.25 
due to the fact that for very hot or very cool stars we did 
not consider their metallicities. Obviously this jump is more 
pronounced for the plots corresponding to [M/H] = —0.7. 

Once the density is known and ag , a\ oggo and a[ M /n] 
have been estimated according to Eq. B5, we find all the 
stars enclosed in each of the eight cubes of size 1.5a P0 . For 
the most populated regions of the parameter space we find 
stars in all these eight cubes. However, if no stars are found 
in a given cube, we simultaneously enlarge its three sizes in 
steps of 0.5x<t po . Since one of the corners of the cube is al- 
ways located on the required parameter point itself, the box 
is enlarged to one side only. In the worst case, this iterative 
procedure ensures finding stars in, at least, three cubes. Fi- 
nally, a similar criterion to that applied for calculating po is 
taken into account for stars cooler than 4000 K and hotter 
than 9000 K, i.e. neglecting their metallicity. 

Having obtained a representative spectrum for each 
cube, S\ 3 , we are now in a position to calculate 
S\(m,t, [M/H]). The weight to be assigned to SV is esti- 
mated as follows 



= 



w 3 



( Pi -<>a \ 2 ( log 8 j -log 8Q A 

V a »a ) e V " losg ° / 



/ [M/H]J-[M/H] \ 2 

V "[M/Hlo J ; (B6) 



Z)j=i a 3 W 3 (logg 3 - log 30) = 
V fc =1 a^([M/H]^-[M/H] )=0, 



(B8) 



where k represents the number of available cubes that con- 
tain stars after the iterative procedure explained above 
(3 < k < 8), and a 3 are coefficients to be applied for cor- 
recting the weights W J to be assigned to the spectra S\ J 
when combining them to obtain S\(m,t, [M/H]). Only val- 
ues in the range < a? < 1 are acceptable for these coef- 
ficients. Among the solutions we choose the one that min- 
imizes the modification of the intrinsic weights W 3 's. This 
is performed by calculating the projection of the solution 
hyperplane onto the point (a 3 = 1, j = 1, k) of the k- 
dimensional space. When no non-negative a 3 solutions can 
be found taking into account the k cubes, we start a proce- 
dure where we neglect the contribution of a given cube in 
Eq. B8. This is performed following the criterion of remov- 
ing first the cubes with lower intrinsic weights according to 
SNjW 3 where SN n j $ is given by 



SN n 



(B9) 



If no satisfactory solutions are found, we select the one 
which minimizes the difference between the required atmo- 
spheric parameters and those resulting from the solution 
star having taken into account a 3 's as in Eq. B8. The final 
spectrum, S\(m,t, [M/H]), is obtained as follows 



S x (m,t,[M/H}) = 



£* ,.s x ; n'ir ; 



(B10) 



where 6 3 , logg^ and [M/H] J are given by 



where S\ 3 are normalized to the selected wavelength ref- 
erence interval AA ro f and scaled according to the derived 
a 3 W 3 . 



APPENDIX C: MODEL UNCERTAINTIES 

In this Appendix we perform a number of tests to evaluate 
the major uncertainties affecting the obtained SSP model 
predictions presented in this paper. We first focus on the 
uncertainties derived from varying the model calculation 
details, particularly in the way that we compute a repre- 
sentative stellar spectrum. We then explore the effects of 
varying the adopted theoretical prescriptions for calculating 
the stellar evolutionary tracks and isochrones. To illustrate 
how significant a given effect is, we have chosen to measure 
the variation of the most important feature in this spectral 
range, i.e. the Call triplet, by means of the CaT* index. 



E - 8 i w i CI Uncertainties derived from the gridding of the 

,JY 3 w j stellar library 

loggj _ Ei=i^ log9 i w j (B7) Here we study the effects of varying the way in which we 

Eiii w i compute a stellar spectrum according to the description 

J^i-! l M / H li W i 



[M/H] 



w 3 



t 



since the most populated cubes would be assigned considerably 
We then solve for the following undetermined linear sys- higher weights, an effect which tends to emphasize the asymmc- 

tem of equations tries present in several regions of the parametrical space. 
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Table CI. Uncertainties of model calculation details on the CaT* index for SSPs of age 12.59 Gyr, unimodal IMF of slope 1.3 and 
different mctallicities 



Test 



Reference Adopted 



CaT* (adopted)-CaT* (reference) (A) 

[M/H]=-0.68 [M/H] = -0.38 [M/H]=0.0 [M/H]=+0.20 



T e ff where [M/H] is neglected 


4000 K 


3750 K 


-0.080 


0.059 


0.046 


0.052 


TpflP where [M/Hl is neglected 


4000 K 


4250 K 


0.020 


-0.005 


-0.012 


0.003 


Y)Pi(Pi < Pm)/J2Pi 


99.7% 


100% 


-0.037 


-0.020 


0.000 


0.006 


22 PiiPi < Pm)/ /J Pi 


99.7% 


90% 


-0.012 


0.015 


0.028 


0.008 


Starting finding box size (xa P0 ) 


1.5 


0.5 


-0.090 


-0.118 


0.000 


0.016 


Starting finding box size (xor po ) 


1.5 


1.0 


-0.071 


-0.020 


0.000 


0.001 


Starting finding box size (x<r P0 ) 


1.5 


3.0 


-0.037 


-0.047 


0.009 


0.033 


Finding box size enlargement per iteration (x<r P0 ) 


0.5 


1.5 


-0.005 


0.002 


0.015 


-0.021 


p range fraction adopted for <r PM 


1/10 


1/6 


-0.166 


0.002 


0.016 


0.059 


p range fraction adopted for <r PM 


1/10 


1/8 


-0.083 


-0.010 


-0.008 


0.030 


p range fraction adopted for cr PM 


1/10 


1/12 


0.000 


0.029 


-0.003 


-0.001 



given in § 3.2 and, particularly, in Appendix B. We se- 
lect for this purpose a representative set of SSP models 
of age 12.59 Gyr (old stellar populations are more sen- 
sitive to the lack of cool stars) and metallicities [M/H] 
= -0.68, -0.38, O.Oand + 0.2. 

We first test to vary the lower temperature limit, below 
which the metallicities of the stars are not taken into account 
when calculating a representative spectrum corresponding 
to a given set of atmospheric parameters. Table CI lists 
the obtained residuals for the CaT* index when increasing 
this temperature from the adopted value, i.e. 4000 K, to 
4250 K, and when decreasing it to 3750 K. We note that 
above 4000 K the obtained differences change very softly. 
Varying this limit towards temperatures lower than 3750 K 
has virtually no effect when compared to the values obtained 
for this temperature limit, due to the fact that most of the 
stars have either solar metallicity or have been assumed to 
be solar when no metallicity values were available in the 
literature. We conclude that this temperature limit is mostly 
affecting the integrated spectrum of SSPs of [M/H] = —0.68, 
but the largest variation in the CaT* index is ~ 1% of the 
CaT* index strength. 

We also explored the effect of varying the density value 
Pm where the cumulative density fraction of stars reaches the 
upper percentile from 99.7 to 100 and 90 (see Table CI). 
In practice for an upper percentile value of 100, this as- 
sumption means that there is a single point in the stellar 
parametric space where a$ m , 0"i ogSm and cr[M/H] m adopt the 
values corresponding to the minimum uncertainties in the 
determination of the stellar parameters. 

Another interesting test is to explore the effect of vary- 
ing the size of the parametrical cubes in which we start to 
find stars. We adopted in this paper 1.5x<t po but we list in 
the same table the values obtained for 0.5 1.0 and 3x<t P q. A 
different, but related test, is to vary the applied enlargement 
to the size of these cubes when no stars are found after each 
iteration. The obtained residuals when adopting 1.5xa P0 , 
rather than 0.5xct po , are tabulated in Table CI. It is clear 
that the starting finding box size plays a more relevant role 
than the way in which the code expands this box when no 



stars are found. The largest residual, i.e. — 0.12A, is found 
for [M/H] = —0.38 when adopting 0.5 x a P0 . However, for 
larger starting size values the obtained differences decrease. 

Finally, we show the effect of adopting for ae M , cr\ oggM 
and o"[m/h] m several fractions of the total range covered by 
these parameters. Our reference value is 1/10, but Table CI 
lists the residuals when we set this limit to 1/6, 1/8 and 
1/12. By increasing this fraction we increase the weights 
assigned to the stars in the less populated regions of the 
parametrical space. This allow us to include stars whose pa- 
rameters are more distant from the requested point. The 
largest fractions (e.g. 1/6) provides the largest effect on the 
Call triplet feature for the smallest and largest metallicities 
(i.e. [M/H] = -0.68 and [M/H] = +0.2), where the CaT* 
weakens by 0.17A for [M/H] = —0.68 and strengthen by 
0.06A for [M/H] = +0.2 with respect to our reference value. 
However, these differences become smaller as we approach 
to our reference value, whilst for smaller fractions (e.g. 1/12) 
the obtained differences are almost negligible for all metal- 
licities. Overall we conclude that the SSP model predictions 
for the largest and smallest metallicities are subject to the 
largest uncertainties. 



C2 Effects of varying the theoretical prescriptions 
of the stellar tracks 

Here we focus on the effects caused by varying the theoret- 
ical prescriptions adopted for calculating the stellar evolu- 
tionary tracks and isochrones. We first discuss the effects 
derived from the most uncertain stellar evolutionary phase, 
i.e. the AGB. Then we discuss the effects of specific theoret- 
ical prescriptions affecting intermediate-aged stellar popula- 
tions, such as the adopted convective overshooting scheme. 
Next we focus on those prescriptions affecting old-aged stel- 
lar populations. Finally we explore the effects of adopting 
stellar isochrones build-up on the basis of non-solar element 
ratios. 

We first show the effects of the AGB stellar evolution- 
ary phase. We illustrate this study by means of performing 
a very crude approach, which is to remove this phase from 
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Table C2. The effect of the AGB on the CaT* for SSPs of dif- 
ferent ages and metallicities. We use models of unimodal IMF of 
slope 1.3. 



Age(Gyr) A C aT * ( AG B [no-ycs] ) (A) 

[M/H]=-0.68 [M/H]=-0.38 [M/H]=0.0 [M/H]=+0.20 

0.20 -0.499 

1.00 -0.087 0.157 

3.16 0.000 0.185 0.120 

10.00 0.000 -0.058 0.035 0.029 



the isochrones. Table C2 lists the obtained residuals for the 
CaT* index for SSPs of different ages and metallicities. For 
old stellar populations the effect of this evolutionary phase 
is very small, and the obtained residuals are of the order of 
the uncertainties derived in Appendix CI. As expected, the 
largest residuals, which can reach as high as ACaT* ~ 0.5 A, 
are obtained for ages below ~ 1 Gyr. These residuals are sig- 
nificantly larger than the ones listed in Table CI. Obviously, 
this is a very crude test, which completely neglects our cur- 
rent knowledge of this late stellar evolutionary phase. In 
practise one expects to obtain smaller residuals due to un- 
certainties in the AGB. This test shows us that the results 
presented in this paper for old stellar populations are stable 
against variations in the input physics of the most uncertain 
stellar evolutionary phase. 

We plot in Figure CI the obtained spectra for an SSP 
of solar metallicity and 0.2 Gyr synthesized on the basis of 
including and excluding the whole AGB stellar evolutionary 
phase. We see that the overall shapes of these spectra are 
remarkably different. For these intermediate-aged SSPs the 
AGB stars are able to shape the continuum around the Ca II 
triplet feature in a way which is characteristic of M-type 
stars. Furthermore, without the AGB phase the Paschen 
lines clearly dominate the synthesized SSP spectra. This ex- 
ercise also shows us to what extent is important to work with 
an appropriate stellar spectral library, with a great coverage 
of the stellar atmospheric parameters including cool stars. 

We have also explored the effects of varying other im- 
portant theoretical prescriptions of the stellar evolutionary 
tracks, such as the adoption of a convective overshooting 
scheme. The new Padova isochrones (GOO) have been up- 
dated on the basis of a significantly milder overshooting than 
that adopted for the old Padova isochrones (B94). Further- 
more, GOO also provide a canonical set of solar metallicity 
isochrones, following the classical Schwarzschild criterion for 
the convective boundaries (i.e. without overshooting). The 
overshooting is only important for stellar masses larger than 
~ 1.2 M©, therefore old-aged stellar populations are insen- 
sitive to the adopted scheme. We have transformed all these 
isochrones to the observational plane following the same pre- 
scriptions described in § 2.3. We then computed SSPs for 
different ages (i.e. from 0.2 to 10 Gyr) for solar metallic- 
ity and unimodal IMF (fx = 1.3). In Table C3 we list the 
obtained CaT* residuals. For SSPs of 0.2 and 1.0 Gyr the 
CaT* values provided by the canonical isochrones are larger 
than those based on the new Padova set, which are larger 
than those obtained from the old Padova isochrones. In fact 
we obtain a total residual ACaT* ~ 1 A. As expected, 
no significant differences are found for old stellar popula- 
tions due to the adopted overshooting scheme (see the first 




8400 8500 8600 8700 8800 8900 

Wavelength (A) 



Figure CI. Synthetic SSP spectra for solar metallicity, 0.2 Gyr 
and unimodal IMF (fi = 1.3). For the upper spectrum we have 
included the AGB phase, whilst for the bottom spectrum we have 
excluded this stellar evolutionary phase. The spectral resolution 
is 1.5 A (FWHM). 



Table C3. CaT* index residuals due to different theoretical pre- 
scriptions in the stellar evolutionary tracks. The obtained resid- 
uals are calculated for SSPs of different ages (indicated in gi- 
gaycars), solar metallicity and unimodal IMF (fi = 1.3). 



Age ACaT*(Overshoot[no-yes]) ACaT* (Padova[old-new] ) 



0.20 


0.481 


-0.600 


1.00 


0.826 


-0.146 


3.16 


-0.016 


-0.631 


10.00 


-0.003 


-0.533 



column of residuals in Table C3), since stars with masses 
smaller than ~ I.IMq present a radiative core during the 
MS. Fig. C2 shows the main differences between these three 
sets of isochrones for different ages. 

For old-aged SSPs we find CaT* residuals larger than 
0.5 A when comparing the old and new Padova isochrones. 
The CaT* values provided by the old Padova set are smaller 
because they predict slightly cooler RGBs as shown in 
Fig. C2. For giant stars with temperatures smaller than 
~ 3500 K the Can triplet decreases with decreasing tem- 
perature (see Paper III). An exhaustive comparison of the 
theoretical parameters used to build-up these two sets of 
isochrones is out of the scope of this paper and we refer the 
interested reader to the GOO paper. However, it is worth re- 
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Figure C2. Plot of the new, old and canonical (i.e. without over- 
shooting scheme) Padova isochrones for solar metallicity and 0.2 
and 10 Gyr. See text for details. 



Table C4. CaT* residuals obtained for SSPs of different ages 
and metallicitics (unimodal IMF of slope 1.3) calculated on the 
basis of the a-enhanced and scaled-solar isochrones of S00. We 
use a unimodal IMF (fi = 1.3). 



Agc(Gyr) 


ACaT*( 


o-enhanced - 


scaled-solar) (A) 


[M/H] = - 


0.38 [M/H] = 


0.0 [M/H] = +0.20 


0.20 




0.954 




1.00 




0.692 


0.555 


3.16 


0.299 


0.475 


0.435 


10.00 


0.269 


0.541 


0.517 



calling that if we replaced the model grids shown in Fig. 15 
with those synthesized on the basis of the old Padova set 
of isochrones, which yield smaller CaT* strengths, we would 
derive smaller IMF slopes, i.e. A[i ~ 1 for the two galaxies 
plotted in this figure (see § 8.2). Therefore, users of these 
SSP model predictions should be warned that, for this rea- 
son, IMF slope estimates must be taken into account on 
a relative basis. Despite the caveats, we believe that our 
predictions are still useful for obtaining robust conclusions, 
particularly when discussing possible trends among galaxies 
(see C03). 

Elliptical galaxies show an enhancement of Mg over 
Fe when compared to scaled-solar SSPs model predictions 



(e.g. Peletier 1989; Worthey et al. 1992; V97; K00; V01A). 
Therefore, it is worth exploring the Call triplet strengths 
in SSPs synthesized on the basis of enhanced a-elements 
mixtures. Such models require both stellar tracks calculated 
with appropriate opacity tables and stellar spectral libraries 
of similar a-enhanced ratios. However, no consensus has 
yet been reached for the degree of enhancement of differ- 
ent a-elements. For example, a constant enhancement for 
each a-element was adopted by VandenBerg et al. (2000) 
and Kim et al. (2002), whilst a somewhat more empirically- 
guided mixtures have been used by SW98 and S00. These 
authors find that, for a given total metallicity, the adoption 
of a-enhanced mixtures leads to higher effective tempera- 
tures with respect to scaled-solar for all evolutionary phases, 
and that, at relatively high metallicities and old ages, an a- 
enhanced isochrone cannot be mimicked by using a scaled- 
solar isochrone of different metallicity. It is worth noting 
that all these stellar models adopt calcium enhanced ratios 
as it is an a-element. However the Ca 4227 line does not 
seem to track Mg but Fe in elliptical galaxies (V97; Trager 
et al. 2000a; TOO; V01A; PS02). On the other hand we do 
not expect any significant variation in the isochrone as a re- 
sult of varying the degree of enhancement of this element, as 
its contribution to the total metallicity is lower than 0.5%. 
The other important problem in synthesizing a-enhanced 
SSPs is the lack of stellar spectral libraries of appropriate a- 
enhanced mixtures. In particular, the empirical library used 
by our model follows the Galactic disk element ratios. 

In spite of the fact that we are unable to build up such 
fully self-consistent a-enhanced SSP models for this spectral 
range, we are in position to evaluate the isochrone effects 
on the synthesized spectra. For this purpose we have made 
use of both a-enhanced and scaled-solar isochrones of S00, 
which are constructed on the basis of almost the same in- 
put physics as that used in GOO. Table C4 lists the residuals 
obtained for the CaT* index for SSPs of different ages and 
total metallicities. We find that the a-enhanced isochrones 
provide larger CaT* strengths for all the tabulated ages and 
metallicities. This is due to the fact that the RGB of these 
isochrones is hotter and therefore there is a smaller number 
of stars falling into the cool temperature regime, where the 
strength of the CaT* decreases as a function of decreasing 
temperature (see § 4.1 and Paper III for more details). We 
also see that the obtained residuals decrease as function of 
increasing age. In particular for old stellar populations of to- 
tal metallicities [M/H] > we obtain ACaT* ~ 0.5 A, whilst 
for [M/H] = -0.38 we obtain ACaT* ~ 0.3 A. Therefore, if 
we plotted these results on Fig. 15 we would obtain steeper 
IMF slopes for these two elliptical galaxies. Put in other way, 
if a Salpeter IMF is used, the a-enhanced isochrones lead to 
a larger calcium underabundance problem. 
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